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Preface 

The  Advanced  Simulation  Concepts  Theory  Program  at  NHL  during  FY  35 
was  focused  toward  applying  the  computational  tools  developed  in  previous 
years  to  provide  guidance  and  supervision  to  the  experiments  on  the  GAMBLE 
II  Facility,  continue  the  work  on  imploding  gas  puff  plasmas,  with  and 
without  a  central  core  plasma,  develop  a  multimaterial  radiation 
hydrodynamic  capability,  initiate  a  program  to  understand  the  non-LTE 
electron  distributions  in  electrical  discharges  as  it  relates  to  and 
impacts  on  the  understanding  and  interpretation  of  experimental  results, 
explore  the  usefulness  of  a  gyrokinetic  corona  model  to  understand  the  flow 
of  current  as  the  plasma  implodes,  to  continue  investigating  the  x-ray 
conversion  efficiency  of  proton  beams  interacting  with  a  planar  target,  and 
finally,  provide  analyses  and  Interpretations  of  the  experimental  data. 

Numerical  simulations  were  performed  to  provide  guidance,  supervision, 
and  interpretlon  for  the  gas  puff  experiments  conducted  on  the  GAMBLE  II 
Facility.  Comparisons  were  made  between  the  implosion  and  radiation 
dynamics  based  on  current  waveforms  representative  of  the  generator  with 
and  without  the  plasma  erosion  opening  switch.  The  numerical  simulations 
were  done  for  both  neon  and  argon  gas  puff  plasmas.  The  predictions  were 
then  borne  out  experimentally  showing  that  neon  was  well  matched  to  the 
GAMBLE  II  generator  In  terms  of  K-shell  x-ray  production  where  Argon  is 
not.  Most  of  the  K-shell  neon  x-ray  emission  resided  in  the  helium-  and 
hydrogen-like  resonance  lines  while  argon  is  a  prolific  L-shell  radiator. 

Investigations  were  also  continued  on  stagnating  an  imploding  gas  puff 
plasma  onto  a  low  density  core  plasma  for  usefulness  as  a  plasma  switch  or 
radiation  enhancement.  It  was  determined  that  the  process  of  stagnation 
has  the  effects  of  stepping-down  the  photon  energy  and  converting  the  K- 
shell  photons  to  L-shell  photons.  The  total  yield  did  not  remain 
constant.  In  fact,  the  presence  of  the  post  consistently  produced  a  higher 
total  yield  in  the  case  of  argon  than  the  corresponding  implosions  in  the 
absence  of  a  central  post. 

A  fully  self-consistent  multi -material  radiation  MKD  model  (DRACHMA) 
has  been  developed  which  includes  detailed  non-LTE  ionization  dynamics  and 
radiation  transport.  An  arbitrary  number  of  materials  can  be  treated 
either  as  mixtures,  discrete  shells,  or  as  a  combihation  of  the  two. 
Numerical  simulations  were  performed  for  an  argon  gas  puff  plasma  imploding 
on  a  neon  core  plasma.  The  numerical  simulations  indicate  that  both 


radiation  enhancement  and  filling  in  of  the  spectra  occurs  under  the  right 
conditions  during  the  stagnation  phase  of  the  implosion. 

A  new  PRS  corona  model  has  been  formulated  on  the  basis  of  gyrokinetic 
theory.  Collisionless  highly  magnetized  plasma  on  the  exterior  of  the  PRS 
load  can  conduct  appreciable  current,  even  when  moving  at  the  E/B  drift 
velocity.  This  gyrokinetic  corona  has  been  examined  in  a  self-similar 
limit  and  the  propagation  theory  developed  into  a  fully  self-consistent 
nonlinear  wave  equation. 

As  an  initial  foray  into  the  domain  of  high  current  opening  switches, 
we  have  developed  the  Eroding  Foil  Switch  (EFS)  model  to  explain  the 
observed  kinematics  of  foils  exposed  to  large  currents  in  a  coaxial 
geometry.  By  calculating  the  mass  erosion  rate  (from  the  upstream  aide  of 
a  foil)  due  to  nonlinear  magnetic  diffusion  and  ohmic  heating,  one  can 
obtain  a  quite  detailed  amd  fairly  accurate  picture  of  the  foil  switch 
conduction  phase.  We  are  involved  in  the  Western  Research  Corp.  (WRC)  foil 
switch  program  as  a  probable  test  bed  of  the  EFS  model. 

The  proton  beam  target  interaction  studies  has  continued.  Numerical 
simulations  were  carried  out  for  beam  intensities  currently  available  with 
existing  pulse  power  generators.  Results  were  obtained  for  a  one  MeV  mono- 
energetic  proton  beam  incident  on  a  planar  iron  target  of  u  8m 
thickness.  The  beam  is  taken  to  be  a  square-shaped  pulse  of  10  nsec 
duration  with  a  flux  of  10^^  protons  cm"^  sec“\  The  simulations  of  the 
target  response  indicate  that  radiation  transport  and  losses  are 
significant  and  influence  the  time  history  of  the  temperature  and  density 
profiles.  Neglect  of  radiation  transport  and  losses  leads,  in  many  cases, 
to  erroneous  behavior  of  the  beam  plasma  interaction. 

A  program  has  been  initiated  to  Investigate  the  non-LTE  electron 
distributions  in  electricad  discharges.  The  electron  energy  distribution 
function  is  the  connecting  link  between  the  electrical  and  radiative 
behavior  of  electrical  discharges;  namely,  it  determines  both  the  strength 
of  the  discharge’s  electrical  and  thermal  conductivities  as  well  as  the 
strength  of  the  ionization  and  excitation  rates,  and  therefore  the 
radiation  loss  rates.  Thus,  it  is  essential  that  electron  distribution 
functions  be  calculated  for  pulse  power  x-ray  machines,  which  should  be 
designed  to  be  as  energy  efficient  as  is  physically  possible  in  converting 
the  heated  free  electron  energy  distributions  to  kilovolt  x-rays.  Without 
a  good  knowledge  of  the  dynamics  of  the  electron  distribution  function,  the 


V 


full  potentialities  of  pulse  power  x-ray  generators  cannot  be  well 
assessed. 


From  the  electrical  and  thermal  conductivities,  one  calculates  the 
current  distributions  (and  thus  the  details  of  the  implosion  dynamics). 
Moreoever,  it  has  long  been  appreciated  that  the  conduction  properties  of 
pulse  power  loads  are  anomalous,  i.e.  they  are  not  understood  in  terms  of 
classically  calculated  Spitzer  conductivities.  Within  the  DOE  IGF 
community,  flux  correctors  have  long  been  in  use  to  sharply  reduce 
classical  heat  conduction  in  steep  temperature  and  density  gradients. 

Experience  has  also  taught  this  community  that  the  values  of  the  flux 
corrector  can  change  markedly  depending  on  experimental  conditions.  Its 
use  also  contributes  to  the  poor  agreement  that  is  usually  obtained  in 
matching  theory  with  detailed  experimental  x-ray  data.  A  discussion  is 

presented  of  the  program  and  its  future  goals. 

Finally,  spectral  analysis  of  GAMBLE  shots  and  comparisons  were  made 
with  previous  BLACKJACK  5  performance.  The  extremely  efficient  performance 
of  GAMBLE-II  (with  respect  to  yield  relative  to  current)  clearly  indicates 
the  possibility  of  upgrading  the  yields  of  the  larger  machines  like 

BLACKJACK  5  without  increasing  their  overall  power.  Employment  of  a  Plasma 

Erosion  Opening  Switch  (PEGS)  on  such  machines  is  one  obviously  promising 
approach.  Current  neon  puff  yields  on  BLACKJACK  5  are  30  KJ ,  about  a 
factor  of  3  higher  than  those  obtained  on  GAMBLE-II  current  of  1 .2  MA, 
Strict  I  scaling  would  suggest  that  hundreds  of  kj  could  be  obtained  with 
neon  puffs  on  BLACKJACK  5.  Clearly,  this  would  violate  conservation  of 
energy  and  such  an  expectation  or  objective  is  therefore  unrealistic. 
However,  the  current  yields  of  30  kj  on  BLACKJACK  5  reflect  an  increase, 
relative  to  GAMBLE-II,  of  less  than  I^. 


In  addition  to  the  work  presented  here  there  have  been  a  number  of 
reports  on  related  activities.  They  include; 

1.  Radiative  Properties  of  Puffed-Gas  Mixtures:  The  Case  of  Optically 

Thick  Plasmas  Composed  of  Two  Elements  with  Similar  Atomic  Numbers", 

J.  P.  Apruzese,  and  J.  Davis,  Journal  of  Applied  Physics  p.  4349, 
Janueiry  1985. 

2.  Opacity  Broadening  as  a  Density  Diagnostics  for  Spot  Spectroscopy",  J. 

Apruzese,  Journal  of  Quantitative  Spectroscopy  and  Radiative  Transfer 
33.  p.  1  February  1985. 

3.  Radiative  Properties  of  Puffer-Gas  Mixtures,  D.  Duston,  J.  Davis  and 

C.  Agritellis,  J.  Appl.  Phys.  57_,  1  Feburary  1985. 

4.  An  Analytic  Voigt  Profile  Escape  Probability  Approximation,  J;  P. 

Apruzese,  J.  of  Quantitative  Spectroscopy  and  Radiative  Transfer, 

3i,  p.  447 

5.  Ion-Beam  Deposition,  Heating  and  Radiation  from  an  Aluminum  Plasma,  J. 

E.  Rogerson,  R.  W.  Clark  and  J.  Davis,  Physicial  Review  A,  Vol.  32, 

May  1985). 

6.  Memo  Report  #5609  "The  Ionization  Equilibrium  of  Optically  Thick  Argon 

Z-Pinch  Plasma  for  Electron  Temperatures  Between  25  and  65  eV,  J.  P. 
Apruzese,  P.  C.  Kepple,  and  J.  Davis,  July  1985. 

7.  Memo  Report  #5590  "Advanced  Concepts  Theory  Annual  Report  1984,  Plasma 

Radiation  Branch,  June  26,  2985. 

3.  Memo  Report  #5613  "Dynamics  of  Imploding  Neon  Gas  Puff  Plasmas  -  I", 

R.  W.  Clark  and  J.  Davis,  July  19,  1985. 

9.  Memo  Report  #5615  "SIMPLODE:  an  Imploding  Gas  Puff  Plasma  Model  I. 

Neon",  J.  Davis,  C.  Agritellis  and  D.  Duston,  July  26,  1985. 

10.  Memo  Report  #  5482  "Hydrodynamic  Response  of  a  Radiation  A1  Plasma 

Excited  by  a  Proton  Seam",  J.  E.  Rogerson,  R.  Clark,  and  J.  Davis, 
December  1984. 


During  FY  85,  GAMBLE-II  at  NRL  was  successfully  operated  both  with  and 
without  the  Plasma  Erosion  Opening  Switch  (PEGS),  using  neon  gas  puff  loads 
created  with  a  2.5  cm  diameter  supersonic  nozzle.  One  task  of  the  Plasma 
Radiation  Branch  related  to  this  work  was  the  analysis  of  the  detailed  neon 
spectra  which  were  emitted  on  these  shots,  thereby  clarifying  the 
differences  in  the  plasmas  produced  with  vs.  without  use  of  the  PEGS.  Vfhen 
these  analyses  were  complete,  a  broad-brush  comparison  of  the  performance 
characteristics  and  efficiency  of  GAMBLE-II  with  that  of  Maxwell 
Laboratories'  BLACKJACK  5  was  undertaken. 


II.  GAMBLE-II  NEON  PUFF  SHOT  ANALYSES 

A  substantial  number  of  successful  neon  puff  shots  occurred  on  GAMBLE- 

II  with  and  without  the  PEGS.  Experimentally,  it  has  been  found  that  use 

of  the  PEGS  eliminates  any  prepulse  and  reduces  the  current  risetime  from 

-  60  ns  to  -  20  ns.  The  price  paid  is  that  typical  peak  currents  with 

the  PEOS  are  0.8  -  0.9  MA,  while  not  utilizing  the  PEGS  results  in  greater 

peak  currents  (-  1.3  MA) .  Given  the  expected  I^  yield  scaling,  one  would 

expect  that  neon  K-shell  yields  with  the  PEOS  would  fall  below  those 

c  0  35  >  4  '' 

obtained  without  the  PEOS  by  a  factor  -  i  ~’~^q j  ••  5  •  However,  yields 
obtained  with  the  PEOS  were  usually  only  a  factor  of  2  below  those  obtained 
without  the  PEOS,  clearly  indicating  greater  overall  efficiency  when  the 
PEOS  is  used.  Analysis  of  prototype  spectra  with  and  without  the  PEOS  was 
undertaken  to  answer  the  question:  what  property  or  properties  of  the 
imploded  plasma  are  responsible  for  the  improvement  of  yield  relative  to 
current  when  the  PEOS  is  utilized? 

Manuscript  approved  February  3, 1986. 
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The  two  Shota  in  which  the  highest  quality  K-shell  spectral  data  were 
obtained  are:  shot  2531  (no  PEGS)  and  shot  2585  (with  PEGS).  They  are 
also  entirely  representative  of  a  much  larger  data  base.  We  therefore 
concentrate  the  spectral  analysis  on  these  shots.  In  each  shot  over  90$  of 
the  radiated  K-ahell  energy  was  contained  in  the  Ly  o  and  He  a  lines.  The 
diameter  in  each  case  was  estimated  as  1.8  mm  from  pinhole  images,  although 
some  overexposure  may  render  this  measurement  only  an  upper  limit.  Knowing 
the  pulse  widths  (29  ns  for  2531  :  21  ns  for  2585),  the  following  measured 
power  outputs  were  obtained. 

TABLE  OF  MEASURED  POWER  OUTPUTS 


Ly  a  POWER 
(Joules/nsec) 


He  a  POWER 
(Joules/nsec) 


PEAK  CURRENT  (MA) 


2531  8.0  26.4  1.30 

(no  PEGS) 

2585  15.0  12.0  O.83 

(with  PEGS) 

Using  the  collisional-radiative-equilibrium  (CRE)  model  of  the  Plasma 
Radiation  Branch,  excellent  spectral  fits  matching  all  the  above  line 
powers  to  better  than  8$  were  found  for  the  plasma  parameters  given 
below.  These  parameters  may  be  interpreted  as  representative  time-  and 
space-averaged  values  for  those  shots. 


TEMPERATURE  (eV) 


ION  DENSITY  (CM"-’) 


2531 

150 

4.0  X  10 

(no  PEOS) 

2585 

135 

3.4  X  10 

(with  PEOS) 

The  higher  temperature  achieved  in  the  PEGS  shot  2585  is  reflected  in 
its  higher  Ly-a  power  as  seen  in  the  first  table  presented  above. 
Therefore,  the  PEGS  -  created  plasmas  seen  to  be  better  neon  K-shell 
radiators  because  they  are  somewhat  hotter  than  non-PEGS  -  created  plasmas 
while  preserving  about  the  same  size  and  density.  The  above  densities 
constitute  nearly  100%  of  the  cold  neon  gas  originally  puffed  in  for  the 
shots.  This  is  in  severe  contrast  to  BLACKJACK  S's  plasmas  analyzed 
previously  by  NHL,  as  seen  below. 


III.  COMPARISON  OF  PREVIOUS  BLACKJACK  5  AND  CURRENT  GAMBLE- II  PERFORMANCE 


Based  on  mature,  detailed  spectroscopic  analyses  of  aluminum/ argon  for 
BLACKJAC  5,  and  neon  presented  above  for  GAMBLE- II,  we  can  now  compare  the 
plasma  parameters  obtained  with  both  machines  as  they  relate  to  the  machine 
parameters. 

The  following  table  is  self-explanatory. 

COMPARISON  OF  BLACKJACK  5  AND  GAMBLE- II 


Machine  or  Plasma 

Property  BLACKJACK  5  GAMBLE- II 


Peak  Current  3-5-5  MA 


0.8  -  1.4  MA 


Current  Risetime 


70  nsec  (main  pulse)  60  ns  (no  PEOS);  20n3  (PEOS) 


Implosion  relative  peak  density 

to  current  risetime  occurs  within  10ns  of 

current  peak 


peak  density  occurs 
1 1 0  nsec  after  current  start 


Loaded  Mass 

fraction  ending  up  20  %  >  70» 


in  :<-3hell  olasma 


Fraction  of  radiation 
which  ia  K-aheLl 


30-50? 


>  70? 


K-Shell  pulse  l0-20n3  20-40n3 

radiating  time 

Ion  densities  few  x  10^^  few  x  10'^ 

achieved 

Temperature  300  eV  -  1 800  eV  100-200  eV 

achieved 

K-3hell  yields  up  to  40  kJ  up  to  4  kJ 

achieved 

It  is  clear  from  the  above  table  that  GAMBLE-II  ultimately 
incorporates  a  much  larger  fraction  of  the  cold  loaded  mass  into  the  K- 
shell  plasma,  and,  correspondingly,  a  much  larger  fraction  of  its  emitted 
radiation  is  X-shell  lines.  Detailed  energy  budgets  confirm  that  a  lot  of 
BLACKJACK  5*3  wires  and  argon  puffs  ended  up  as  L-shell  ions  emitting  L- 
sheli  radiation. 

IV.  IMPLICATIONS  FOR  DNA  OBJECTIVES 

The  extremely  efficient  performance  of  GAMBLE-II  (with  respect  to 
yield  relative  to  current)  clearly  indicates  th?  possibility  of  upgrading 
the  yields  of  the  larger  machines  like  BLACKJACK  5  without  increasing  t.heir 
overall  power.  Employment  of  a  Plasma  Erosion  Opening  Switch  (PEGS)  on 
such  machines  is  one  obviously  promising  approach.  Current  neon  puff 
yields  on  BLACKJACK  5  are  30  KJ,  about  a  factor  of  3  higher  than  those 
obtained  on  GAMBLE-II.  However,  with  a  peak  BLACKJACK  current  of  3.7  MA, 
as  opposed  to  GAMBLE-II  current  of  1.2  MA,  strict  I"*  scaling  would  suggest 
that  hundreds  of  kJ  could  be  obtained  with  neon  puffs  on  BLACKJACK  5. 
Clearly,  this  would  violate  conservation  of  energy  and  sucn  an  expectation 
or  objective  is  therefore  unrealistic.  However,  the  current  yields  of  30 
kJ  on  BLACKJACK  5  reflect  an  increase,  relative  to  GAMBLE-II,  of  less  than 


I  .  Certainly  these  yields  can  be  improved,  and  we  strongly  reconnnend 
systematic  experimentation  with  different  nozzle  radii  and  plenum 
pressures. 
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B.  SIMPLODE;  A  Non-LTE  Dynamic  Pinch  Model  of 

Imploding  Gas  Puff  Plasmas.  II.  Argon 

I.  INTRODUCTION 

In  an  earlier  paper^  results  were  presented  on  the  radiation  dynamics 
of  an  imploding  neon  gas  puff  plasma  for  conditions  typical  of  the  GAMBLE 
II  pulse  power  generator  at  NRL.  It  was  shown  that  for  a  variety  of 
initial  conditions  the  neon  gas  puff  is  a  well  matched  load  for  GAMBLE  II 
producing  upwards  of  6  kilojoules  of  K-shell  radiation,  with  the  bulk  of 
the  line  radiation  emanating  from  the  resonance  lines  of  hydrogen-  and 
helium-like  neon.  It  was  further  observed  that  the  peak  of  the  K-shell 
radiative  yield  occurred  well  after  the  peak  in  the  current  pulse 
suggesting  that  it  may  be  possible  to  enhance  the  radiative  yield  by 
optimizing  the  implosion  dynamics.  Although  the  yields  are  strongly 
dependent  on  peak  current,  results  indicate  that  application  of  the  PEGS 
(Plasma  Erosion  Opening  Switch)  improves  the  quality  of  the  implosion, 
i.e.,  the  re-assembled  plasma  exhibits  less  flaming  and  appears  more 
uniform.  However,  further  experimentation  with  the  PEGS  is  necessary 
before  we  can  say  with  certainty  that  it  further  enhances  the  K-shell  yield 
and  improves  significantly  the  quality  of  the  implosion. 

In  this  paper  we  focus  our  attention  on  argon  gas  puff  plasmas. 
Interest  in  argon  is  twofold:  (1 )  provide  estimates  of  the  radiative  yield 
for  a  variety  of  puff  gas  and  machine  conditions  and  (2)  investigate  the 
feasibility  of  creating  conditions  favorable  to  supporting  a  population 
inversion  and  lasing  in  the  3i~3i'  transitions  in  neon-like  argon. 
Preliminary  calculations  indicate  that  it  may  be  possible  to  generate 
conditions  favorable  to  lasing.  These  results  are  being  reported  on 
elsewhere."  The  objective  of  this  investigation  is  to  characterize  the 
radiation  dynamics  of  an  imploding  argon  gas  puff  plasma  in  much  the  same 
fashion  as  we  did  for  neon.  Unlike  neon,  argon  is  not  a  particularly  well 
matched  load  to  the  GAMBLE  II  generator  for  producing  prolific  K-shell 
yields.  However,  the  disappearance  of  K-shell  radiation  will  more  than  be 
made  up  for  by  the  appearance  of  prolific  L-shell  yields.  All  the 
calculations  are  made  with  the  SIMPLODE  code  containing  an  argon  non-LTE 
ionization  dynamic  model  which  has  been  reported  on  elsewhere." 


vie; 
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II.  RESULTS  AND  DISCUSSION 


For  reference  purposes  we  will  present  results  showing  how  some  of  the 

key  plasma  parameters  vary  in  time  as  the  implosion  evolves.  The 

standarized  case  is  represented  by  a  30  ugm/cm  argon  gas  puff  plasma  with 
initial  outer  (inner)  radius  of  1.55  cm  (0.95  cm)  and  a  total  column  length 
of  4  cm.  The  mass  distribution  is  uniform  and  homogeneous.  The  current 
waveform  driving  the  implosion  is  shown  in  Figure  1  and  has  a  peak  value  of 
0.875  megamps.  The  shape  of  the  current  waveform  is  typical  of  the  GAMBLE 
II  generator  with  the  PEGS  with  the  exception  of  the  tail  beyond  about  150 
nsec  which  has  been  installed  for  computational  expediency.  The 

experimentally  measured  current  may  have  a  longer  decay  time.  However, 
since  it  is  the  peak  current  that  dictates  much  of  the  implosion  dynamics 
of  interest  here,  we  shall  ignore  the  differences.  In  Figure  2  the 

temporal  history  of  the  inner  (A)  and  outer  (B)  radius  is  shown.  Note  that 
the  inner  radius  collapses  on  axis  in  about  125  nsec  and  that  the  outer 
radius  pinches  to  0.18  cm  (which  is  about  1/10  the  initial  radius)  and 
exhibits  a  bounce  at  about  200  nsec  indicating  that  both  temperature  and 
density  achieve  values  high  enough  to  produce  a  gas  pressure  that  retards 
and  exceeds  the  magnetic  pressure.  Correspondingly,  the  velocity  reaches  a 
maximum  value  of  13cm/usec  at  150  nsec  and  becomes  zero  at  the  bounce 
point,  as  shown  in  Figure  3.  The  ion  density  is  shown  in  Figure  4  and 
indicates  a  peak  ion  density  of  about  5  x  10^^  cm~^  at  the  time  of  peak 
compression  -  the  bounce  point.  Similarly,  in  Figure  5  the  temporal 
behavior  of  the  temperature  is  presented.  The  temperature  peaks  at  about 
170  nsec  with  a  maximum  value  of  about  180  ev.  The  total  radiative  yield 
as  a  function  of  time  over  the  4  cm  column  length  is  shown  in  Figure  6. 
The  peak  value  is  about  1.7  kjoules  and  occurs  roughly  at  the  bounce 
time.  Note  that  the  temperature  and  density  peak  at  different  times.  This 
is  caused  by  radiative  cooling  as  the  plasma  runs  in.  This  can  be  seen  by 
comparing  with  the  radiative  yield  curve  in  Figure  5.  The  plasma  starts 
radiating  significantly  at  about  130  nsec  and  has  a  total  radiation  pulse 
of  about  40  nsec.  The  overwhelming  contribution  .o  the  radiative  yield 
comes  from  the  L-shell  and  will  be  shown  and  discussed  in  detail  later. 
These  Figures  summarize  the  behavior  of  some  of  the  '<ey  plasma  parameters 
characterizing  the  implosion  and  radiative  dynamics  associated  with  our 
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reference  case.  It  is  apparent  from  this  reference  case,  which  will  be 
shown  to  be  a  judicious  choice,  that  argon  is  not  well  suited  to  GAMBLE  II 
in  terms  of  producing  significant  K-shell  radiation  but  is  well  suited  for 
generating  interesting  levels  of  L-shell  radiation. 

We  will  now  summarize  the  results  of  a  series  of  numerical  simulations 
fixing  either  the  mass  or  the  thickness  of  the  plasma  sheath, 
Ar  •  Rg  -  R^.  The  simulations  were  performed  for  two  typical  current 
waveforms:  the  no-switch  case,  which  is  the  current  waveform  associated 

with  the  GAMBLE  II  generator  without  the  PEGS,  having  a  peak  current  of 
1.25  Megamp  which  is  achieved  in  70  ns  and  is  shown  in  Figure  7,  and  the 
current  waveform  with  the  PEGS  shown  in  Figure  i  and  already  discussed  in 
the  context  of  Figures  2-6.  Because  of  the  voluminous  amount  of  data  it 
was  decided  that  the  best  way  to  present  the  results  was  to  choose  the 
maximum  value  of  the  appropriate  variable  from  each  computer  run  and  plot 
it  as  a  Tunction  of  mass  per  unit  length  or  sheath  thickness.  For  example, 
for  fixed  Ar,  corresponding  to  Rg  -  R^  (•  1.55  “  0.95)  •  0.60,  the  maximum 
velocity,  is  shown  in  Figure  8  as  a  function  of  mass  per  unit  length 

for  both  current  waveforms.  In  comparing  the  differences  between  the  two 
current  cases  it  is  evident  that  the  higher  peak  current  case  gives  rise  to 
higher  peak  implosion  velocities.  The  higher  the  current  the  greater  the 
magnetic  acceleration  the  faster  the  implosion  velocity  for  a  fixed  mass  of 
plasma.  In  Figures  9  and  10  are  displayed  the  peak  temperatures  and  ion 
densities.  For  low  values  of  M/I  the  kinetic  energy  of  runin  is 
sufficiently  large  such  that  as  the  plasma  reassembles  it  thermalizes  at  a 
high  temperature  which  causes  the  following  events:  (1)  the  plasma  over 
heats  causing  a  burn  through  of  the  ionization  stages  and  radiates  as  a 
bremsstrahlung  emitter,  which  is  not  very  efficient  in  cooling  the  plasma 
and  (2)  because  of  the  high  temperature  the  back  pressure  is  high  at 
relatively  low  densities  impeding  the  implosion  and  produces  to  a  "hollow 
implosion",  i.e.,  the  plasma  cannot  compress  to  relatively  high  densities 
which  accounts  for  the  behavior  shown  in  Figures  9  and  10.  As  M/l 
increases  the  driving  current  becomes  less  effective  in  generating  high 
implosion  velocities  which  means  there  is  less  energy  available  for 
conversion  to  thermal  energy.  Hence,  the  plasma  temperature  is  reduced  and 
the  implosion  proceeds  as  shown  in  the  Figures  9  and  '0.  The  peak  density 
exhibits  broad  profile  reaching  a  value  of  approximately  '0  ^  ions/cm-’  in 


the  switch  case.  The  behavior  of  the  peak  ion  density  as  a  function  of 
M/i  can  be  explained  in  terms  of  the  driving  current.  In  the  switch  case 
as  M/Z  increases,  the  peak  density  first  increases,  then  levels  off  and 
subsequently  decreases.  This  behavior  can  be  understood  as  follows;  for 
low  and  moderate  values  of  M/1  the  current  is  over-  and  well  -  matched  to 
the  load,  respectively,  in  terms  of  the  velocity,  temperature,  and  achieved 
density.  However,  as  M/1  increases  the  mismatch  between  current  and  load 
also  increases  producing  slower  implosion  velocities  and  subsequently 
reduced  peak  densities.  In  the  no-switch  case  very  large  values  of  M/1 
were  not  considered,  so  the  effect  is  not  shown.  However,  for  large  enough 
M/1  similar  behavior  to  the  switch  case  would  be  observed.  The  peak  K- 
and  L-shell  radiative  yield  as  a  function  of  M/1  is  shown  in  Figure  11  for 
Ar  -  0.60.  Since  the  yield  is  a  strong  function  of  density  and  tempera¬ 
ture  the  bdiavior  of  the  curves  is  self-evident.  The  peak  yield  with  the 
switch  achieves  a  maximum  value  of  about  3.8  kJ  at  approximately  22  ygm/cm 
while  the  no  switch  case  has  a  maximum  value  of  about  6  kJ  at  about 
35  ygm/cm  for  a  cm  long  plasma  column.  These  values  are  in  reasonably 
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good  agreement  with  the  experimental  observations. 

Another  series  of  numerical  simulations  were  performed  varying  Ar  for 
fixed  M/1.  The  outer  radius  was  chosen  fixed  at  1.55  cm  while  the  inner 
radius  was  permitted  to  vary  from  0.55  o  1.25  cm.  Tne  fixed  value  for  M/1 
was  chosen  as  30  ugm/cm.  Figures  12-15  show  the  peak  values  of  velocity, 
temperature,  density,  and  yield,  respectively,  as  a  function  of  Ar.  The 
peak  velocity  decreases  as  Ar  increases  for  both  current  waveforms.  This 
simply  reflects  the  fact  that  for  small  Ar  the  inner  boundary  has  further 
to  go  and  experiences  little  or  no  back  pressure  whereas  for  large  Ar  the 
inner  boundary  reaches  the  axis  early  in  time  and  allows  buildup  of  back 
pressure  which  inhibits  the  implosion  velocity.  Hence,  the  peak  velocity 
is  greater  for  smaller  Ar.  Peak  temperatures  and  densities  are  shown  in 
Figures  13  and  11*  respectively.  The  peak  temperature  in  the  no  switch  case 
decreases  as  Ar  increases  mainly  because  the  quality  of  the  implosion 
deteriorates  resulting  in  a  cooler  plasma.  This  behavior  is  tied  to  the 
lack  of  a  good  compression  as  shown  in  Figure  1U  for  the  peak  ion  density 
as  Ar  increases  for  fixed  mass  and  jxB  force.  In  the  switch  case  the 
peak  temperature  remains  relatively  constant  over  the  range  of  Ar  values. 
This  just  reflects  the  fact  that  it  is  possible  to  heat  a  portion  of  the 
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pla3ma  to  about  180  eV  -  the  plasma  doesn't  get  sufficiently  hot  enough  to 
radiate  prolifioally  so  that  the  temperature  is  not  significantly  reduced 
by  radiation  cooling.  This  is  further  borne  out  in  Figure  1^  where  the 
peak  ion  density  decreases  as  Ar  increases,  indicating  a  reduction  in  the 
peak  ion  density,  thereby  reducing  the  number  of  emitters  available  to  cool 
the  plasma.  In  Figure  the  peak  yields  for  this  series  of  simulations  is 

presented.  As  Ar  increases  the  no  switch  peak  yield  increases,  peaks,  and 
then  decreases  whereas  the  switch  peak  yield  steadily  decreases  as 
increases  indicating  the  the  best  peak  yields  for  this  case  occur  when  the 
temperature  and  density  are  also  high.  Note,  however,  that  the  yield  scale 
(vertical  axis)  does  not  vary  greatly  ~  in  fact  the  difference  between  the 
highest  and  lowest  values  differ  by  about  10X  for  both  current  waveforms 
indicating  that  the  initial  implosion  conditions  conspired  in  such  a  way  as 
to  generate  peak  values  of  temperature  and  density  so  as  to  produce  little 
variation  in  the  peak  yield. 

Figure  16  displays  the  temporal  variation  of  the  L-  and  K-shell 
radiated  power  for  the  M/J.  -  30  ugm/cm  and  Rg(a^)  »  1  .55  (0.95)  cm  case 
superimposed  on  the  PEOS  current  waveform.  The  peaks  in  the  thermal  L-and 
K-shell  emission  profiles  occur  at  slightly  different  times  with  the  K- 
shell  peak  just  preceding  the  L-shell  peak.  This  is  in  qualitative 
agreement  with  the  observations.^  However,  the  K-shell  radiative  yield  is 
signif iceintly  greater  in  the  experiments  than  predicted  by  SIMPLODE.  A 
possible  explanation  for  this  disparity  is  that  the  experimental  K-shell 
yield  is  due,  for  the  most  part,  to  nonthermal  energetic  electrons  as 
suggested  by  Jones  and  Kania.^  The  L-shell  accounts  for  most  of  the 
radiated  power  and  peaks  around  2  terrawatts.  The  K-shell  radiated  power 
for  these  conditions  is  in  the  noise.  Mote  also  that  the  peak  radiated 
power  comes  well  after  the  peak  in  the  current  pulse.  This  behavior  was 
also  observed  in  the  neon  implosions  for  current  waveform  having  sharp 
risetimes.  Apparently  the  plasma  absorbs  the  energy  and  just  percolates 
for  a  time  and  then  coasts  in  until  it  reassembles  on  axis. 

Finally,  in  Figure  17  the  total  emission  spectra  is  shown  for 
M/I  *  35  ygm/cm  and  R  (Rq),  -  1  .55  (0.95)  cm  for  the  switch  case  at  i77 

O  n 

nsec  into  the  run,  which  is  near  the  peak  in  temperature.  The  ion  density 
at  this  time  is  about  l  .Sx'O  ^  cm~^  and  the  average  charge  state  is 
'2.56.  The  final  pinched  radius  is  0.28  cm.  All  of  the  lines  appearing  in 
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the  spectra  are  L'Shell  lines  with  the  peak  emission  occurring  at  bands 
centered  around  20,  70,  and  MOO  eV.  The  free-bound  continuum  around  1  keV 
falls  off  rapidly  and  the  K-shell  thermal  emission  is  virtually  extinct  at 
the  temperature  associated  with  the  implosion. 

Ill .  SUMMARY 

A  number  of  numerical  simulations  using  the  SIMPLODE  code  were 

performed  to  describe  the  behavior  of  an  argon  gas  puff  plasma  using 
machine  parameters  typical  of  the  NHL  GAMBLE  II  pulse  power  facility.  The 
SIMPLODE  code  is  a  0-D  non-LTE  dynamic  pinch  model  and  should  be 

interpreted  as  a  "scopping"  code  and  is  not  presented  as  gospel.  Its 

virtues  are  a  valid  treatment  of  the  non-LTE  ionization  dynamics  and 
radiation,  phenmonologlcal  representation  of  the  implosion  dynamics,  and 
rapid  running  on  a  Vax  computer.  Its  principal  deficiency  is  that  some  of 
the  relevant  physics  is  deleted  such  as  field  and  current  penetration, 
development  and  growth  of  Instabilities,  and  plasma  transport  properties, 
to  mention  a  few.  However,  it  must  be  said  that  in  comparison  with  more 
sophisticated  multidimensional  radiation  MHD  codes  such  as  PRISIM,  SIMPLODE 
simulates  the  overall  morphology  of  the  implosion  physics  with  the 

exception  of  instabilities  and  provides  reasonable  agreement  with  the 
radiation  yields.  It  is  in  this  context  that  SIMPLODE  is  applied. 

The  results  of  the  simulations  shows  that  argon  gais  puff  plasmas  can 
be  heated  and  compressed  to  temperatures  and  densities,  respectively,  with 
the  GAMBLE  II  generator,  that  produce  a  good  plasma  radiation  source  in  the 
soft  x-ray  region  of  the  spectra  with  yields  of  the  order  of  2-2.5 
terrawatts  over  a  30-40  ns  pulse.  For  fixed  dr  -  0.60,  the  yield  has  a 
peak  value  of  about  6  kJ  for  an  M/I  -  32  ugm/cm  for  the  no  switch  case  and 
a  peak  yield  of  about  4  kJ  for  an  M/I  »  20  ugn/cm  for  the  switch  case. 
For  a  fixed  M/I  -  30,  the  peak  yield  as  a  function  of  dr  between  0.30  and 
0.90  cm  maintaining  the  outer  radius  fixed  at  1.55  cm  remains  in  the 
neighborhood  of  5  to  6  kJ  for  the  no  switch  case  and  between  '  .5  to  i  .75  '<J 
for  the  switch  case.  Obviously  higher  yields  favor  lower  M/l  in  the 
switch  case  because  of  the  lower  peak  current. 

In  conclusion,  even  with  the  limited  parameter  space  investigated, 
imploding  an  argon  gas  puff  plasma  on  the  GAMBLE  II  generator  produces  an 
efficient  L-shell  soft  x-ray  plasma  radiator  source. 
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.  PRISM;  Plasma  Radiation  Implosion  Source  Model 

I,  STAGNATING  ARGON  PLASMAS 

One  of  the  problems  of  interest  is  that  of  a  stagnating  plasma  in 
which  a  puff  gas  implodes  onto  a  dense  core  or  post  plasma.  These 
configurations  allow  for  a  relatively  uniform  implosion  along  the  axis.  A 
stagnating  plasma  also  allows  one  to  control  in  a  crude  manner  the  yield 
characteristics  of  the  implosion.  That  is,  if  one  compares  implosions  in 
which  a  post  is  present  and  implosions  in  which  only  a  background  gas  is 
present,  acre  L-shell  radiation  occurs  when  a  post  is  present.  Similarly, 
provided  the  implosion  velocity  is  low  enough,  more  K-shell  radiation 
occurs  when  there  is  no  post.  Although,  the  goal  of  the  PRS  program  is  to 
maximize  the  harder  K-shell  radiation,  some  applications  may  require 
outputs  in  the  L-shell.  The  concept  of  a  post  thus  allows  one  more  degree 
of  flexibility  in  matching  the  machine  and  mass  loading  to  the  desired 
spectral  output. 

In  the  process  of  stagnation,  the  kinetic  energy  of  the  inflowing 
material  is  converted  into  internal  energy.  In  turn,  the  internal  energy 
is  partitioned  between  ionization  and  thermal  energy.  Some  of  this  energy 
is  radiated  away  from  the  system  or  is  carried  to  another  portion  of  the 
plasma.  A  shock  wave  develops  and  propagates  in  the  direction  opposite  to 
the  inflowing  material  and  defines  the  boundary  between  the  region  of  hot 
dense  stagnated  matter  from  the  region  of  the  relatively  cool,  and  less 
dense  imploding  matter.  In  the  idealized  case  (Y  -  5/3),  the  ra.io  of  the 
densities  is  16  for  a  cylindrical  implosion.  If  a  core  plasma  is  present, 
some  of  the  kinetic  energy  is  taken  up  by  the  core  which  subsequently 

implodes  and  bounces  off  the  center  of  the  column  at  r-0.  I^epending  on  the 

conditions  present,  the  core  can  interact  with  the  imploding  material  after 
bouncing  off  the  center  axis  before  all  of  this  material  has  become 
shocked. 

Results  for  argon  and  neon  stagnating  plasmas  have  been  obtained.  The 
previous  report  presented  the  neon  results.^  In  this  report,  we  will 
briefly  describe  the  argon  results.  The  design  which  we  have  modeled 
contains  some  important  simplifications.  First,  both  the  core  (or  ''post") 

and  the  imploding  puff  are  taken  to  be  a  single  material.  Second,  an 
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idealized  initialization  is  used  in  that  a  low  density  puff  is  initially 
given  a  uniform  velocity  and  is  allowed  to  drive  into  the  post  plasma  which 
is  at  a  much  higher  density.  This  is  in  lieu  of  a  driving  magnetic 
field.  This  approximation  allows  us  to  isolate  the  implosion  physics  from 
the  physics  of  the  particular  machine  which  might  be  under  consideration. 

As  in  our  previous  work,^  we  have  performed  simulations  with  and 
without  central  posts.  The  post  consists  of  an  initially  cold  (0,02  eV) 
argon  gas  at  a  density  of  5  x  gm/cm^  and  a  radius  of  1  mm.  The  argon 

puff  gas  which  implodes  on  the  core  has  an  initial  density  of  5  x 

gra/cm^  and  extends  out  to  1.5  cm.  Tne  puff  gas  starts  with  a  'uniform 

velocity  across  its  width.  Figure  1  shows  an  implosion  with  a  core  present 

and  a  puff  gas  velocity  of  U  x  10^  cm/sec.  Four  times  are  shown:  4  nsec, 
10  nsec,  20  nsec,  and  60  nsec.  At  4  nsec,  contact  has  been  made  and  the 
leading  edge  of  the  puff  gas  has  been  heated  to  about  800  eV,  At  10  nsec, 
the  stagnating  shock  can  be  seen  propagating  outward  with  a  density  plateau 
developing  at  the  interface.  A  shock  can  also  be  seen  moving  into  the  core 
plasma.  Note  that  the  scale  for  the  radius  is  given  in  log-iQ  cm  so  that 

the  core  can  be  distinguished.  By  20  nsec,  the  shock  in  the  core  is  well 

developed  while  the  outward  propagating  shock  has  almost  emerged  from  the 

puff.  At  60  nsec,  the  entire  configuration  is  in  the  process  of 

disassembling,  thought  not  shown  here,  the  shock  in  the  core  plasma 
rebounds  off  the  center  axis  at  about  30  nsec. 

If  we  now  examine  the  radiative  yields  for  cases  with  and  without 

central  core  plasmas  at  several  velocities  there  are  some  interesting 
features  which  did  not  appear  in  our  earlier  work  with  neon.  Because  of 
the  higher  Z  of  argon  we  exended  the  velocity  range  out  to  8  x  10' 

cm/sec.  Our  previous  work  on  neon  (Z  »  10  as  opposed  to  Z  «  18  for  argon) 

had  only  covered  the  range  from  1  x  10  cm/sec  to  3  x  10'  cm/ sec.  The 

first  feature  which  is  different  is  t.he  divergence  of  the  total  yield  as 
velocity  increases  as  shown  in  Fig.  2.  At  all  velocities,  the  cases  with  a 
core  plasma  present  (solid  circles)  showed  a  greater  yield.  The  difference 
between  the  two  cases  is  very  much  evident  at  3  x  10^  cm/sec.  The  line 
yield  for  the  L-shell  is  shown  in  Fig.  3.  Again,  the  divergence  appears 
and  in  fact  Fig.  3  very  closely  mimics  Fig.  2.  The  strong  turnover  in  the 
open  points  (no  core  plasma  present),  is  due  to  burning  through  the  L-shell 
ionization  stages  at  velocities  greater  than  x  iQ'  cm/sec.  The  K-line 
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yield  ia  shown  in  Fig.  4.  At  the  highest  velocity,  8  x  10'  cm/sec,  the  K- 
llne  yield  has  flattened  out  when  a  core  is  not  present  while  it  has 
increased  dramatically  when  a  core  is  present.  This  reflects  the  fact  that 
with  the  absence  of  a  core  more  of  the  kinetic  energy  of  the  puff  gas  is 
converted  to  thermal  energy  with  the  result  that  even  the  He-like 
ionization  stage  is  being  burned  through. 

The  results  for  argon  at  these  higher  velocities  show  that  there  can 
be  substantial  differences  in  the  total  radiative  output  when  a  post  is 
present.  This  is  a  different  result  than  that  which  was  obtained  from  neon 
implosions  initiated  with  velocities  in  the  1  -  3  x  lo”^  cm/sec  range.  At 
puff  gas  velocities  greater  than  4  x  To”^  cm/sec,  the  yields  with  a  post 
present  are  significantly  larger  for  the  K-shell  yield  as  well  as  for  the 
total  and  L-shell  yields. 
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Figure  '  -  Density  and  temperature  profiles  for  an  implosion  in  which  a 
core  is  present  and  an  initial  puff  gas  velocity  of  4  x  ’0'  cm/sec.  The 
densities  are  multiplied  by  a  factor  of  10^  gm/cm^  and  the  temperatures  are 
in  eV.  Four  times  are  shown:  4  nsec,  10  nsec,  20  nsec,  and  60  nsec. 
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VELOCITY  (  10 '  CM/SEC  ) 

Figure  2  -  Total  radiation  yield  for  the  argon  inipiosions  as  a  function  of 
initial  velocity.  Closed  circles  represent  cases  where  a  core  is  present 
and  open  circles  represent  cases  where  there  is  no  core. 


D .  Multimaterial  Radiation-Hydrodynamic  Modelinj 


A  fully  aelf-conaistent  multimaterial  radiation  hydrodynamic  model  has 
been  developed  for  application  to  problems  of  interest  to  DNA.  It  features 
detailed  non-LTE  ionization  dynamics  and  radiation  transport  using  a 
probabilistic  formalism.  An  arbitrary  number  of  materials  can  be  treated 
in  full  generality  -  as  a  mixture,  as  discrete  shells,  or  as  a  combination 
of  these.  A  multi -material  puff -gas  implosion  calculation  will  be 
described,  involving  an  argon  puff -gas  incident  upon  a  neon  core  plasma. 

I .  INTRODUCTION 

A  number  of  recent  DNA  pulsed-power  experiments  have  made  use  of 
mixtures  or  layered  configurations  consisting  of  multiple  materials.  For 
example,  experiments  have  been  performed  involving  annular  puff  gasses 
imploding  onto  dissimilar  axial  puff  gasses  or  low  density  plasma  cores. 
In  addition,  a  number  of  Interesting  x-ray  laser  schemes  propose  to  use 
photons  from  a  strong  transition  in  one  element  to  create  a  population 
inversion  in  another.  There  are  many  important  applications  where  a 
general  multimaterial  radiation  hydrodynamic  algorithm  is  necessary.  For 
this  reason,  the  DRACHMA  code  (detailed-configuration  radiation-transport 
and  coupled  hydrodynamics  multimaterial -mixture  algorithm)  has  been 
developed. 

As  a  test  of  the  multimaterial  algorithm,  and  to  demonstrate  its 
capabilities,  a  simulation  of  an  annular  argon  puff-gas  imploding  onto  a 
central  neon  core  plasma  was  performed.  Sandia  National  Laboratory  has 
carried  out  a  program  involving  puff -gas  Implosions  onto  low  density  foams 
to  determine  the  feasibility  of  creating  a  homogeneous,  uniform  gain  medium 
along  the  central  core  axis.  Preliminary  experimental  results  indicate 
high  radiative  conversion  efficiencies  from  relatively  uniform  reproducible 
implosions.  The  geometry  of  the  simulation  corresponds  to  the  Sandia 
experiments,  but  the  materials  chosen  do  not.  They  were  chosen  for 
convenience  and  for  ease  of  comparison  with  existing  single  material 
calculations. 


II.  THEORETICAL  MODEL 


Since  the  implosion  of  a  cylindrical  puff-gas  plasma  can  result  in  a 
substantial  fraction  of  the  total  plasma  energy  being  radiated  away,  the 
hydrodynamic  evolution  of  the  plasma  may  be  appreciably  altered.  The 
interaction  of  the  plasma  dynamics,  the  atomic  physics  of  the  plasma 
components  and  the  transport  of  radiation  is  both  nonlinear  and  nonlocal 
and  must  be  calculated  self -consistently. 

Discussion  of  the  theoretical  coodel  will  be  separated  for  convenience 
as  follows:  (a)  hydrodynamics  and  thermal  conduction  (b)  ionization  and 
atomic  physics  and  (c)  radiation  emission  and  transport. 

(a)  Hydrodynamics  and  Thermal  Conduction 

The  basic  hydrodynamic  variables  of  total  mass,  momentum,  and  total 
energy  are  transported  in  one  dimension  using  a  numerical  scheme  with  a 
sliding-zone  version  of  flux-corrected  transport.  A  special  gridding 
algorithm  is  employed  which  moves  zones  in  a  Lagrangian  fashion  and  adjusts 
the  mesh  in  order  to  resolve  steep  gradients  in  the  flow.  The  hydrodynamic 
equations  solved  are 
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where  p  is  the  mass  density  summed  over  species, 

P  -  I  N  m^  (U) 

3  ®  ® 

With  Ng  the  number  density  of  species  s  and  mg  the  associated  mass,  u  is 
velocity,  P  is  pressure,  E  is  total  energy  density,  E  .  is  the  rate  of 


energy  loss  or  gain  due  to  radiation,  and  K  is  the  thermal  conductivity. 
The  thermal  conduction  is  calculated  implicitly,  using  an  iterative  Crank- 
Nicholson  scheme. 


Since  density  did  not  exceed  solid  density  in  this  problem,  a  simple 
equation  of  state  was  assumed, 

P  -  2/3  (E  -  1/2  pu^  -  E^^^)  (5) 

where  is  the  potential  energy  due  to  ionization  and  excitation,  summed 

over  species.  A  single  temperature  model  was  employed. 

kT  -  P  /  ♦  I  H)  (6) 


where  T  is  temperature.  and  Ng,  the  electron  number  density,  are 

calculated  from  the  ionization-radiation  equations  which  are  described 
below.  A  single  temperature  approximation  is  valid  in  the  core  plasma, 
where  the  equilibrium  time  is  of  the  order  of  picoseconds,  and  it  is 
adequate  in  the  stagnation  region,  where  the  equilibration  time  can  be  of 
the  order  of  nanoseconds.  In  the  blowoff  plasma,  it  is  a  marginal 
approximation,  but  the  consequences  are  minor  since  little  radiation  is 
emitted  from  this  region,  and  most  of  the  energy  is  carried  by  the 
electrons  in  the  blowoff. 

The  local  rate  of  change  of  energy  due  to  radiation  transport,  E^^^  » 
will  be  discussed  below. 


(b)  Ionization  and  Atomic  Physics 


The  ionic  populations  in  the  plasma  are  determined  by  a  set  of  atomic 
rate  equations  of  the  form 


df . 

d't  ’  ^  '^Jis  ^J3  ~  ^  ”ij3  ^is 


(7) 


i3  is  the  fractional  population  of  atomic  level  i  for  the  plasma 
species  s.  is  the  net  reaction  rate  describing  the  transition  from 


initial  state  j  to  final  state  i.  An  equation  of  this  type  is  constructed 
for  each  of  the  levels  included  in  the  model  for  each  of  the  atomic 
species. 

For  sufficiently  dense  plasmas,  the  effective  populating  and 
depopulating  rates  are  sufficiently  fast  compared  with  the  hydrodynamic 
timescales  that  an  equilibrium  assumption  can  be  justified,  which  involves 
dropping  the  explicit  time  dependence  in  Eq.  (7).  The  plasma  is  then  said 
to  be  in  collisional-radiative  equilibrium  (CRE),  whereby  the  plasma 
ionization  state  responds  instantaneously  to  changes  in  hydrodynamic 
q'uantities. 

Ti.a  rate  coefficients  which  are  used  to  calculate  the  populating  and 
depopulating  rates,  are  calculated  using  various  atomic  calculational 

J  *3 

methods.  The  processes  included  and  the  methods  employed  in  generating  the 
corresponding  rate  coefficients  are  summarized  elsewhere. 

Once  the  set  of  rate  equations,  including  the  optical  pumping  from  the 
radiation  field,  has  been  solved  for  the  level  populations,  f^g,  the 
electron  density  can  be  calculated, 

“e  ■  J  “3  i  '■13  ^3  '3> 
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where  Zj^g  is  the  ionic  charge  of  level  i  for  species  s,  and  Ng  is  the  total 
ion  density  of  species  s. 

The  ionization  and  excitation  energy  can  also  be  calculated  by 
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where  is  the  energy  of  level  i,  measured  from  the  ground  state  of  the 
neutral  atom  of  species  s. 

For  the  simulation  described  below,  the  atomic  model  for  neon  contains 


27  atomic  levels  and  13  emissions  lines,  and  the  model  for  argon  contains 
3^  atomic  levels  and  emission  lines. 


(c)  Radiation  Emission  and  Transport 


Radiation  emission  from  and  absorption  by  a  plasma  are  dependent  on 
the  local  atomic  level  population  densities.  Except  for  optically  thin 

plasmas,  however,  the  level  populations  depend  on  the  radiation  field, 
since  optical  pumping  via  photo  ionization  and  photoexcitation  can  produce 
significant  population  redistribution.  Thus,  the  ionization  and  radiative 
processes  are  strongly  coupled  and  must  be  solved  self -consistently.  In 
this  model,  an  Iterative  procedure  is  used,  where  level  populations  for 
each  material  are  calculated  using  the  radiation  field  from  the  previous 
iteration,  then  using  these  populations  to  calculate  a  new  radiation  field 
and  recalculating  populations  until  convergence  is  achieved. 

A  probabilistic  radiation  transport  scheme  is  employed,  which  forms 
local  angle  and  frequency  averaged  escape  probabilities  for  each  emission 
line  and  for  each  bound-free  process.  This  amounts  to  characterizing  each 
emission  line  and  each  bound-free  process  by  a  frequency  at  which  the 

probability  of  escape  is  maximized  for  the  greatest  number  of  photons, 
since  a  decrease  in  opacity  away  from  the  line  center  or  the  edge  frequency 
can  be  accompanied  by  a  corresponding  decrease  in  emissivity.  Free-free 
radiation  is  treated  with  a  multi frequency  transport  formalism.  The 
radiation  transport  and  emission  spectra  are  calculated  from  these  escape 
probabilities.  The  method  can  treat  comprehensive  atomic  models  and 
provides  good  overall  energetics,  but  cannot  calculate  accurately  certain 
spectral  details  and  lines  with  very  high  optical  depths. 

Inner-shell  opacities  are  included  in  the  model,  since  these  processes 
are  very  importauit  in  cool,  dense  plasma  regions.  Inner-shell  photo¬ 

ionization  cross  sections  for  each  of  the  neutral  elements  are  taken  from 
the  fits  of  Biggs  and  Lighthill,  and  t.he  positions  of  the  ionization- 

dependent  absorption  edges  are  taken  from  the  Hartree-Fock  calculations  of 
Clement i  and  Roetti. 

The  local  rate  of  change  of  energy  in  zone  j,  due  to  radiation 
transport  is  given  by 

S  --II(F  -Tc  F 
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where  Fpi^g  is  the  rate  of  energy  loss  in  zone  k  due  to  a  discrete  radiative 
process  (or  frequency  group)  p  for  species  s,  and  Cpi^jg  is  the  radiative 
coupling  of  zone  k  to  zone  J  for  that  process.  The  couplings  are  functions 
of  opacity,  integrated  over  process  and  photon  path.  In  the  probabilistic 
model,  a  matrix  of  couplings  must  be  computed  for  each  bound-bound,  bound- 
free  and  free-free  process.  In  this  way,  the  net  cooling  and  heating  by 
radiation  emission  and  absorption  between  the  various  aones  of  the  plasma 
is  accurately  taken  into  account. 

III.  RESULTS 

A  simulation  was  performed  for  a  cylindrical  annular  argon  puff-gas  of 
density  5  x  10”®  g/cm,  with  inner  radius  0.55  cm  and  outer  radius  1 .95  cm, 
imploding  radially  at  a  velocity  of  5  x  lo"^  cm/sec  onto  a  central  neon  core 
plasma  of  density  5  x  10”^  g/cm  and  radius  0.10  cm.  The  initial 
configuration  is  shown  in  Fig.  1.  The  temperature  of  the  puff-gas  was 
taken  to  be  about  2.5  eV  initially,  and  a  tenuous  background  argon  plasma 
was  placed  between  the  puff -gas  and  core  plasmas  of  density  5x10'  g/cm. 

These  initial  conditions  were  chosen  to  correspond  approximately  to 
simulations  conducted  with  a  neon  puff-gas  stagnating  on  a  neon  core 
plasma,  described  in  an  earlier  report;  however,  the  implosion  velocity 
has  been  increased  in  order  to  produce  some  K-shell  argon  radiation. 

For  the  first  10  nanoseconds,  the  argon  puff -gas  essentially  coasts 
radially  inward.  The  forward  edge  of  the  puff-gas  is  heated  to  about  100 
eV  through  accretion  of  the  background  plasma.  At  about  10  nanoseconds 
(Fig.  2),  the  puff-gas  makes  full  contact  with  the  core,  and  conversion  of 
kinetic  energy  to  thermal  energy  begins  to  build  an  overpressure  at  the 
interface. 

I  3y  19.4  nanoseconds  (Fig.  3)  the  peak  temperature  near  the  interface 

is  approaching  a  keV,  and  a  narrow  annulus  of  intense  net  radiative 
emission  centered  at  the  contact  region  is  formed.  The  emission  spectrum 
at  this  time  is  shown  in  Fig.  4.  It  is  dominated  by  neon  and  argon  L-shell 
I  radiation  at  this  time,  and  the  emission  falls  off  rapidly  above  1  keV. 

j  The  individual  contributions  of  the  argon  and  neon  components  to  the  total 

I  emission  spectrum  can  be  seen  in  Figs.  5  and  6,  respectively.  In  Fig.  5, 

.  the  argon  spectrum  is  plotted  together  with  the  total  free-free 
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radiation.  The  total  argon  emission  is  substantially  larger  than  the  neon 

emission,  shown  in  Fig.  6,  since  the  bulk  of  the  core  plasma  is  still  very 

cold.  Optical  depth  through  the  plasma  is  plotted  as  a  function  of  photon 
energy  in  Fig.  7.  At  19. nanoseconds,  inner-shell  opacity  is  the  major 
component,  and  since  most  of  the  plasma  is  still  relatively  cool,  inner- 
shell  optical  depths  have  not  dropped  substantially  from  cold-target  values 
except  at  low  energies.  The  neon  cold-target  L-shell  and  K-shell 

absorption  edges  occur  at  0.018  keV  and  0.867  keV,  whereas  the 
corresponding  argon  L-  and  K-shell  edges  occur  at  0.2U5  and  3.203  keV, 
respectively. 

3y  19.4  nanoseconds,  the  overpressure  at  the  interface  has  caused 
shocks  to  propagate  radially  inward  in  the  core  and  outward  in  the  puff- 
gas.  Heating  from  the  soft  x-rays  produced  in  the  interface  region  is 
evident  in  the  dense  core,  and  radiation  heating  is  warming  the  puff-gas. 
Thermal  conduction  is  ineffective  in  transporting  energy  into  the  core 
plasma. 

The  inward  propagating  shock  reaches  the  origin  at  about  3*  nano¬ 

seconds.  Plasma  profiles  at  this  time  are  shown  in  Fig.  8.  Although  the 
neon  plasma  approaches  solid  density  on  axis,  the  temperature  there  is 
still  only  about  50  eV.  This  region  of  very  dense  plasma  radiates  strongly 
in  the  L-shell,  but  its  volume  is  small.  A  broader  annulus  of  hot  plasma 
situated  outside  the  core  plasma  is  producing  most  of  the  radiation.  This 
can  be  seen  from  the  emission  spectrum  at  this  time,  shown  in  Fig.  9,  where 
argon  L-shell  radiation  still  dominates. 

Optical  depths  are  plotted  as  a  function  of  photon  energy  for  this 
time  in  Fig.  10.  The  neon  inner-shell  optical  depths  are  characteristic  of 
a  strong  temperature  gradient,  resulting  in  an  array  of  ionization 
states.  The  ionization-shifted  L-shell  absorption  edges  of  neon  are  seen 
as  a  series  of  steps  below  about  200  eV.  The  increase  in  optical  depth 
above  the  shifted  L-shell  edge  of  neon  is  due  to  the  compression  of  the 
core. 


At  60  nanoseconds  (Fig.  11),  the  entire  plasma  is  moving  radially 
outwards.  The  temperature  of  the  neon  core  plasma  has  fallen  to  about  20 
eV,  and  the  argon  puff-gas  has  cooled  to  about  120  eV. 

The  energy  history  of  the  plasma  throughout  the  first  HO  nanoseconds 
is  shown  in  Fig.  12.  Initially,  about  98%  of  the  total  energy  is  kinetic; 
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the  remainder  is  in  thermal  and  ionization  energies.  The  energy  lost  to 
radiation  begins  to  increase  sharply  at  about  12  nanoseconds,  as  the  argon 
puff-gas  comes  into  contact  with  the  neon  core.  Thermal  and  ionization 
energies  also  begin  to  climb  at  this  time.  Kinetic  energy  goes  through  a 
well  defined  minimum  at  about  38  nanoseconds,  at  which  time  more  than  half 
of  the  total  plasma  energy  has  been  radiated  away.  By  110  nanoseconds, 
more  than  70?  of  the  total  energy  has  been  converted  to  radiation. 

IV.  CONCLUSIONS 

Puff-gas  implosions  show  promise  as  a  means  of  e'‘ficiently  converting 
kinetic  energy  to  radiation.  With  a  central  stagnation  core  plasma,  higher 
plasma  densities  and  total  radiative  yields  can  be  achieved. 

This  simulation  served  as  a  testbed  for  the  DRACHMA  multi-mater ial 
code,  it  accurately  calculated  the  generation  and  transport  of  radiation  in 
a  multi-species  environment.  Including  a  detailed  accounting  of  continuum 
and  inner-shell  opacities. 

Since  radiation  plays  such  a  crucial  role,  it  must  be  generated  and 
transported  accurately.  As  we  have  shown  previously,  optically  thin 
radiation  ignores  the  importance  of  transport  of  radiative  energy  and 
grossly  overestimates  losses.  LTE  or  corona  equilibrium  models  for  the 
level  populations  will  also  lead  to  inaccurate  energetics.  Finally, 
because  the  radiative  energy  loss  and  the  transport  of  energy  in  the  plasma 
substantially  modify  the  hydrodynamic  response,  a  self-consistent 
calculation,  linking  radiation  and  fluid  transport  with  the  ionization 
dynamics,  must  be  performed. 


E.  GYROKINETIC  PLASMA  RADIATION  SOURCE 
CORONA  MODEL 


1.  Corona  Plasma  Phenomenology 

Becanse  pure  magnet ohydro dynamic  models  generally  predict  faster, 
denser  Z-pinch  implosions  than  one  observes  in  pnlse  power  diodes,  a 
common  corrective  hypothesis  has  been  one  form  or  another  of  "anoma¬ 
lous  resistivity* .  In  early  NRL  modela^'^  of  the  plasma  radiation 
source  (PBS),  an  exterior  plasma  doiaain  characterized  by  thin  density 
and  high  temperature  was  invoked  as  a  mechanism  for  softening  the  pinch 
by  providing  an  alternate  current  path.  By  carrying  a  drift  speed  lim¬ 
ited  current  density  and  being  too  tenuous  to  lose  heat  through  ra¬ 
diation,  such  a  plasma  could  perhaps  absorb  a  modest  fraction  of  the 
diode  current  and  thus  rob  the  interior,  denser  load  plasma  of  its 
driving  force,  the  "magnetic  piston".  If  so  the  limitations  of  "pure* 
MED  models  could  be  understood  from  a  more  fundamental  viewpoint. 

Such  a  simple  model  for  "anomalous  resistivity*  proved  severely  lim¬ 
ited.  In  its  original  form,  the  constraints  of  quasi-static  pressure 
balance,  rather  than  smooth  radial  acceleration,  and  a  conventional 
Ohm’s  law,  produced  a  corona  of  consistently  small  current  and  mass, 
relative  to  the  main  load.  On  the  other  hand,  it  is  now  apparent  that 
a  more  general  picture  of  the  exterior  plasma  does  provide  a  poten¬ 
tially  significant  current  channel.  The  basis  for  this  current  lies 
in  two  phenomena  —  the  heating  of  tenuous  gas  by  the  first  part  of 
the  electromagnetic  driving  pulse,  and  the  magnetization  of  the  re¬ 
sulting  plasma  by  the  interior  current  driving  the  main  load  mass. 


A. 


Electron  "Driftavay* 


The  detaile  of  early  heating  are  at  present  unexplored  in  the  PRS 
coamunity.  Such  thin,  hot  plaeaae  ae  aight  occur  in  the  exterior  re¬ 
gion  are  difficult  to  diagnoee  apectroacopically  due  to  the  largo  flu- 
ence  of  XUV  and  eoft  Z-raye  eaitted  by  the  denser  portione  of  the  load 
as  it  begins  to  absorb  the  electroaagnetic  pulse.  Tet  it  is  this  very 
fact  which  aakee  reasonable  the  partial  or  perhaps  eoaplste  photoion- 
ization  of  the  tenuous  load  gas  or  "pnap  oil"  surrounding  the  load. 

The  tendency  of  aagnetically  insulated  transaission  line  (MITL)  sys- 
teas  to  propagate  poorly  at  background  prasenros  anch  exceeding 
10~'*  torr  substaatiatas  this  conjecture^  .  One  can  generally  expect 
atoaic  densitiao  in  the  range  of  10^^  to  10^^  par  ca~’  in  the  ia- 
aediata  exterior  neighborhood  of  a  PRS  load.  When  this  gas  is  par¬ 
tially  ionized  its  radial  flow  velocity  is  initially  saall.  any  large 
axial  E  field  such  as  that  of  the  incoaing  driver  pulse  will  force 
further  ohaic  heating,  and  possibly  drive  unstable  aicroturbulance, 
which  will  aaka  the  medina  even  aore  resistive.  This  heating  will 
proceed  until  the  plasaa  can  no  longer  absorb  energy  through  collisions. 
Hence  a  reasonably  credible  production  scenario*  for  diffuse  exterior 
plasaa  is: 

i.  Early  ohaic  heating  of  tenuous,  partially  photoionized  load  gas 
or  "puap  oil"  drives  the  electron  collision  time  above  the  PRS  ia- 
plosion  timescale  of  [10  -  100]  ns. 

it.  Non-maxwellian  distributions  of  electrons  and  some  lighter  ions 
decouple  from  the  colder  population  and  effectively  drop  the  con¬ 
ductivity,  further  increasing  the  heating  rata. 


Hi.  Haatiag  procaadi  natil  an  aaargy  danalty  ia  raachad  for  which 
all  tha  alactroaa  ara  lockad  np  in  tha  non-maxwallian  population. 

(V.  High  anargy  alactrona  ara  eontinnonaly  captnrad  in  ExB  drift 
aotiona  rathar  than  axial  runaway. 

Tha  dataila  ara  aeeaaaibla  through  a  Fokkar-Planck  calculation,  which 
ia  in  davalopaant. 

B.  Tanuoua  Plaaaa  Corona  —  Gyrokinatic  Curranta  and  Drifta 

Onca  tha  colliaionlaaa  populationa  ara  "boilad  off  of  tha  background 
gaa,  tha  riaing  intarior  currant  carriad  by  tha  dioda  load  tanda  to 
pravant  runaway  ia  tha  uaual  aanaa.  Tha  magnetic  field  demanda  an  avar 
more  iaaulatad  atata  and  aome  fraction  of  tha  tanuoua  gaa  will  drift 
inward  aa  a  magnetinad  plaama.  Tha  magnetic  inanlation  of  thia 
ExB  plaama  ia  not  perfect,  however,  bacanaa  work  muat  be  dona  on 
it  to  make  tha  radial  drift  follow  tha  boundary  conditiona  impoaad  by 
tha  pulaa  power. 

The  power  abaorbed  manifaata  itaalf  in  current  danaity  and  E  fiald 
attenuation,  the  "corona"  therefore  preaenta  an  energy  aink  between 
the  pulaa  power  aourca  and  tha  PRS  load.  Such  magnetized  drift  cur¬ 
rent  danaitiea  have  three  componenta  whan  deacribed  by  a  "gyrokinetic" 
theory.  Tha  firat  ariaea  in  radially  accelerating  tha  guiding  cen- 
tara  due  to  incraaaaa  in  E/B  during  the  driver  pulaa.  Tha  aecond, 
from  tha  azimuthal  angular  momentum  aaaociatad  with  guiding 
cantar  circulation  around  tha  dioda  axia.  Tha  third,  from  the  "gy¬ 
ration”  angular  momentum  about  tha  guiding  canter,  uaually  axpreaaed 


in  tarmi  of  a  magnotic  momant.  la  moving  tha  magnatizad  plasma  to¬ 
ward  tha  dioda  axis  tliaas  coasarvad  angular  momanta  forca  an  incraaas 
in  tlia  partiela  apaad  witli  dacraasing  radius  or  incraasiag  magnatic 
fiald.  Tba  work  raquirad  to  incraasa  this  spaad  comas  from  the  an- 
argy  in  tha  alactromagnetie  drivar. 

Tka  laval  of  tkis  powar  loss  doaa  not  diminish  as  tha  plasma  com- 
prassas  (in  fact  the  magnetization  criteria  improve  near  tha  origin) . 
Moreover,  tha  angular  momentum  conservation  in  tha  tenuous  region  would 
require  an  infinite  amount  of  energy  to  gat  the  magnatizad  plasma  to 
the  origin.  Latar  in  time,  only  tha  onset  of  collisionality, . or  the 
relaxation  of  tha  magnatic  fiald,  can  destroy  tha  constants  of  the  mo¬ 
tion  and  permit  a  dissipative  and  to  tha  process.  Tha  theory  of  gy- 
rokinatic  coronal  phenomena  is  thus  dominated  by  three  considerations: 

t.  Collisionlaaa,  magnatizad  electrons  and  ions  persist  in  gy- 
rokinetic  drifts,  cf.  Fig.l  n-T  contours  of  (lira  ,  tha  number 
of  gyrations  in  a  collision  time; 

it.  Unmagnetized  ions  follow  tha  magnatizad,  gyrokinatic  species 
through  space  charge  fields;  any  regions  of  charge  separation  al¬ 
low  narrow,  intansa  ExB  currant  channels; 

tit.  Lata-tima  populations  may  reach  the  core  plasma,  releasing  stored 
energy,  or  providing  exterior  currant  channels  which  soften  tha  core 
implosion. 


Fig.  1  Contours  of  Qiru  for  at  1.5  10^  Gauss 


ti.  Unmagnetized  ions  follow  the  magnetized,  gyrokinetic  species  through  space  charge 
fields;  any  regions  of  charge  separation  allow  narrow,  intense  E  x  B  current  channels; 

Hi.  Late-time  populations  may  reach  the  core  plasma,  releasing  stored  energy,  or  providing 
exterior  current  channels  which  soften  the  core  implosion. 


II.  Gyrokinetic  Plasma  Model 

The  aatnral  substrata  on  vbich  to  build  a  corona  nodal  is  collision- 
lass  gyrokinatic  theory  and,  as  in  any  Vlasov  description,  tha  fun- 
danantal  class  of  functions  is  that  of  nataaquilibria.  Thasa  tine  sta 
tionary  solutions  of  tha  gyrokinatic  Vlasov  equation  are  associated 
with  tha  free  streaming  particla  orbits  in  the  self-consistent  diode 
alsctroaagnatic  fields.  In  the  following  development  are:  (i)  tha 
nature  of  such  nataaquilibria,  (ii)  tha  minimum  set  of  observable  mo- 
mants  required  to  pose  a  drift  currant  constituitiva  relation,  and 
(iii)  a  general  progagation  aquation  appropriate  to  one  dimensional 
absorption  and  reflection  of  transverse  magnetic  modes.  Tha  later  sec 
tions  show  a  particular  solution  of  the  propagation  problem  and  dis¬ 
cuss  tha  implications  for  tha  PRS. 

A.  Geometry  and  Fundamentals 

Specialiming  the  gyrokinatic  picture  of  Bernstein  and  Catte’  to  a 
double  disk  configuration  with  an  interior  axial  current,  a  general 
calculation  of  electromagnetic  power  flow  across  an  exterior,  colli¬ 
sionless  plasma  annulus  begins  with  the  observation  that  any  ambipo- 
lar  electric  field  present  can  be  represented  by  a  potential  function 
long  as  its  gradients  are  orthogonal  to  the  purely  az¬ 
imuthal  magnetic  field.  This  potential  function  can  be  determined  in 
principle  by  momentum  conservation  between  the  magnetized  species  and 
the  unmagnetized  ions,  if  any.  In  the  simplest  form  of  the  problem 
all  species  are  assumed  magnetized  and  the  radial  field  can  be 
neglected  or  assumed  to  come  to  a  quasi-equilibrium  configuration  just 
strong  enough  to  preserve  the  local  charge  balance.  In  such  a  case 

the  primitive  variables  of  the  model  are  easy  to  summarize. 


la  tha  traasTaraa  magaatic  fiald:  B  =  B9[r,t)  9  ,  tha  TEM  aoda 
alactric  fiald:  E2(r,  f)  =  £{5,  r,  f)  -  Er.4m6(-,  r,  f)  =  E  +  .  rap- 

rasaata  tba  racaaa  aolatioa  of  tha  radial  liaa.  Thia  aoda  aat  iapliaa 
a  radial  drift  valocity:  U  =  — r(Ei/Bj)c  ,  aad  a  total  drift  va- 
loeit7:  V  =  U  +  vhieh  iaeladaa  eircalatioa  about  tha  axia. 

la  additioa.  tha  paeuliar  valocity:  =  Wji  +  (wr  -  ?7r)r  ,  ia  or- 

thogoaal  to  tha  local  gyratioa  radiaa  vaetor.  which  ia  ia  tura  locatad 
at  a  gyrophaaa  aagla:  ^  —  {7,)/Ui)  .  Tha  total  particla 

valocity  ia  thaa  v  as  V  +  . 

Tha  goal  of  tha  gyrokiaatic  thaory  ia  tha  dataraiaatioa  of  oaa  body 
diatributioa  fuactioaa  which  daacriba  tha  colliaioalaaa  aieroatataa 
of  tha  dioda  coroaa.  Siaca  tha  coroaa  productioa  procaaa  ia  aot  ax- 
pactad  to  favor  aay  gyrophaaa,  tha  gyrophaaa  ia  avaragad  ia  dafiaiag 
tha  "gyrokiaatic  gyroavaraga*  diatributioa  fuactioaa, 

/<»,a(r,/i<,,«,0  =  (l/23r)  j  dtp  /a(r, /ia,  «, 

whara  tha  aora  accurata  gyrophaaa  dapaadaat  variablaa  (  /ia,u  )  ara 
aaad,  at  laaat  foraally.  Thia  iapliaa,  howavar,  that  aoma  Vlaaov 
aataaquilibria,  furthar  apacialiaatioaa  of  tha  ,  will  dapaad  ia 

affact  oa  avaragad  valuaa  of  tha  highar  ordar  iavariaata  which  caa  ba 
labalad  by  tha  spatially  local  valuaa  of  lowaat  ordar  iavariaata. 

ia  a  furthar  aiaplificatiou  to  tha  gyrophaaa  avaragad  kiaatic 

equatioa  which  provldaa  the  aataaquilibria,  oaa  haa,  for  thia  situ- 
atioa,  the  additioaal  raault  ^ • (V  x  V)  =  0  . 


All  &ppropria.t«  aataaqnilibria  bara  ara  fanetiona  oi  strict  eoastaata 
of  tha  motion,  aneb  na 


Pa, a  =  mo,  r  u  (1  - 


Wrt  eoa{  ^  ) 
fl  r 


) 


mo,  r  , 


axpraaaad  in  tarma  of  both  gyropbnaa  dapandant  and  indapandant  vari- 
ablaa.  Altarnataly,  an  ineoaplataly  kaovn  eonatant  of  tba  motion  ad- 
mita  bigbar  ordar  approximationa  in  tba  form  of  gyropbaaa  dapandant 
invarianta,  a.g. 


iia 


ntg 

2Be 


{i-6) 


Yitb  a  =  (e,  t)  and 


S  =  2  ^  coa{  0  )  +  (f^)  »»n(  2^  ) 

Wrg  2  ll  o$ 


Vban  tba  "gyraklnatie  gyroavaraga*  oparation  ia  appliad  to  any  func¬ 
tion  of  tbaaa  strict  or  approximate  invarianta,  tba  strict  invariants 
sucb  aa  tba  angular  momantum,  ,  remain  nsaful  argumants;  but  tba 

gyropbaaa  dependant  invariants  ara  affactivaly  reduced  to  tba  ordi¬ 
nary  ones  —  magnetic  momant,  Ha  —  initial  guid¬ 

ing  canter  position,  .  ITota  that  tbia  reduction  in  complexity  is 
not  a  radnetion  in  precision  bacauaa  tba  lower  order  variables  are  re¬ 
ally  just  labels  for  gyroavaragad  bigbar  ordar  invarianta. 


All  sucb  mutually  indapandant  invariants,  numbers  labeling  guiding 
cantar  orbits,  tbarafors  datarmina  (tbrougb  ordar  l/H  )  gyrokinatic 
Vlasov  aetaaquilibria  of  tba  ganaral  form 


fo,a{r,v,  t)  =  ( 


m: 


2t 


}  Po,a[/^a  ?  Pa, a  '  ''o) 


Tba  factor  involving  r  ,  ma  .  and  B9  arisaa  from  tha  velocity  space 
transformation  to  )Ua,  and  ■  Various  functions  Po,a  are  poa- 

sibla,  but  only  tba  lowest  velocity  moments  ara  naadad  to  datarmina 
tba  gyrokinatic  danaity  evolution.  Sucb  momants  ara  time  independent 
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•xpactmtioa  valuta  of  higli  ordar  gyrophasa  avaragad  iavarianta,  ob- 
aarvatioaa  or  calcnlationa  of  tba  corona  prodnctioa  phasa  must  astab- 
liah  tham. 

Tbara  ia  also  no  aagaatic  diffnaioa  ia  tha  gyrokinatic  plasma,  viz. 

atB=:V  X  (U  X  B), 

forcaa  tha  aagaatic  fiald  avolntioa  to  coaaarva  tha  aagaatic  flnx 
tbrongb  any  cloaad  cnrva  that  follows  tba  drift  velocity  Ur  aad  tbe 

power  absorption  is  therefor a  dataraiaad  by  the  dynamics  of  tha  fiald 

E.  . 

B.  Drift  Currant  Density 

Tha  enrraat  daasity  for  tha  gyrokiaatic  fluid  is  not  accessible 
through  a  linear  Ohm's  law.  Tha  axial  drifts  which  constitute  such  a 
curraat  are  iabalancaa  la  aaaa  charged  particle  flow  caused  by  chaagas 
ia  tha  local  macroscopic  fields  aad  tha  spatial  transport  of  guiding 
canters.  Vhaa  tha  gyrokinatic  aquation  is  transformed  to  a  saquanca 
of  moment  equations,  this  curraat  can  be  datarminad  through  the  con** 
straints  imposed  by  momentum  conservation,  and  the  time  invariant  moments 


ifl=  ^  f  di.,dP,,„ 

''  P  J  ml 

and 

p  J  ma 

summed  on  a  tha  index  of  magnetized  species.  The  mass  density 

P  “  mg,  not  “  oIq,  [  (  )  f  dfi(jidP$,a  fo,a  ]  > 

mi  r  J 


ia  of  cooroo  dotorminod  by  partielo  eeasorvatioa.  Tbooo  aomonto  aloao 
dofiao  tho  divorgoaco  of  tbo  prootaro  toasor.  yioldiag  a  carroat  doa- 
■ity 

=  C  ^'‘§7*  *57'  i  ^  I’ 

aad  aa  oqnivaloat  momentum  conservation  equation 

fi 

p  Dt  Ur  =  — jj-  -  r~^dr{r  p  B9  W^)  -  Jz  B9  . 

For  tboso  ioa  apoeiaa  wbicb  arc  aot  aagaotizod,  tbo  aoaoataa.  coaaor* 
▼atioa  rulo  bocoffloa  aodifiad  by  tbo  additioa  of  aa  aabipolar  atraaa 
tara,  a  aealar  proaanra  tara,  aad  tba  aao  of  ,  tha  "eaatar  of 

aaaa"  voloeity,  ratbar  tbaa  aa  tba  flow  fiald,  viz. 

P  Dt  Wg  =  —  <r  drXamb  "*  ^r(B)  —  Jion,z  B$  C  ^  . 

Wo  ia  tbaa  dataraiaad  priaarily  by  tba  aabipolar  atraoaoa  ariaiag 
vbaa  ligbtar  apaeiaa  drift  iaward  abaad  of  tba  baaviar  ioaa.  It  will 
taad  to  aatcb  tba  drift,  with  aoaa  tiaa  dalay,  aa  tba  baaviar 

apaeiaa  ara  alao  implodad. 

C.  Gaaaral  Propagatioa  Eqaatioa 

Siaca  tba  earraat  daaaity  caa  ba  aaaaaaad  only  from  tba  dyaasica  of 
Ur  and  fo,a  •  tba  moat  ganaral  formnlatioa  ia  to  aliminata  J;  from 
tba  field  aqnationa  and  develop  a  wave  aquation  for  Es  to  completa 
tba  picture.  Tba  aaaa  daaaity  and  magnetic  fiald  are  convaniantly  da- 
acribad  by  a  Lagrangian  poaitioa  E(ro,  t)  for  tba  guiding  canter,  vbila 
tba  fiald  aquation  ia  naturally  poaad  in  taraa  of  dimenaionlaaa  vari- 
ablaa  r  =  tf to  ,  r  =  rflg  ,  lo  —  c  to  .  For  example,  one  might  aet 
to  aa  tba  current  riaatima,  providing  a  uaeful  acala  for  Eg  aa 
^z,o  —  ^diode  ^peakl^o  • 


Tlia  Alfvaa  ipaad  ca  =  ariaat  as  a  natural ,  variable 

wave  speed  and  using  tbs  Jacobian  associated  with  tbe  Lagrangian  flow, 
dR{foif‘)Mfo  t  the  position,  density  and  magnetic  field  are  given  by 


r  - 

R{fo,r)  =  fo  -  J  dTx  (i2(ro,ri),ri)  , 


A  aK 


and 


aK 


With  the  r  notation  suppressed  and  the  fields  understood  as  scaled 
by  Eg^o  •  the  material  derivative  is  expanded  and  combined  with  the 
"frozen*  magnetic  field  constraint,  the  current  density  =  4r<o  Jz 
now  becomes 


h  =  ^\9rE.  - 


Wv 


BbW, 


whara  1-*  ,  aad  '  •  A  nonlinaar  wava  aqua- 

tiea  eaa  ba  davalopad  for  Eg  by  aliaiaatiag  tba  aagaatic  fiald  in 
Faraday's  law  and  collaetiag  tboaa  taraa  in  drJz  which  dapand  on  tha 
diaplaeaaaat  currant.  Tha  propagation  aquation  is  than 


r-‘a,(r9,£.)  -  (1  +  -  a,{^)  d,E.  +  d,\^d,&\  = 

c\  c\  e.  Be 


dMrpB,W^)\ 


This  propagation  aquation  points  out  thraa  iaportant  charactaris- 
tica  of  tha  gyrokinatic  population  — 


t.  Tha  Eg  fiold  ovolution  ranains  hyperbolic  in  charactar,  but  tha 
alactromagnatic  phasa  valocity  slows  to  tha  Alfvsn  spaad,  viz.  Vph  — 
(iwVc»)t/a  •  Alfvan  spaad  equals  tha  light  spaad,  tha 

phasa  valocity  is  still  half  tha  vacuum  valua:  tha  Alfvan  spaad  must 
ba  much  largor  than  tha  light  spaad  in  order  to  racover  tha  usual 
macroscopic  transit  times.  This  provides  a  strong  modification  to 
tha  magnetic  diffusion  picture,  where  prompt  wave  propagation  is  usu¬ 
ally  invoked  to  neglect  displacement  current.  Tha  gyrokinatic  medium 
shows  a  natural  pulse  compression  mechanism  which  may  find  useful 
applications . 

it.  Comparing  the  coefficient  of  tha  first  order  time  derivative  of 
Eg  with  the  mass  density  and  magnetic  fiald  evolution  rules,  one 
notaa  that  damping  occurs  when  tha  plasma  compression  exceeds  tha 
magnetic  field  energy  compression,  since  drlnic^ /c'\)  =  -drln{^^)  . 


This  is  just  tlis  conditioa  for  whieli  tlis  guiding  cantors  nra  boiag 
scealarstad  by  sa  iacrosaiag  Eg  fiald.  In  contrast,  a  plasma  vhicb 
is  dacalarating  its  guiding  cantors  will  ganarata  Eg  fiald  through 
this  nachanisa,  a  kind  of  nagativa  conductivity.  Iha  anorgy  flow 
botvaan  particlaa  and  Eg  fiald  is  smoothly  raciprocal  and  non- 
dissipativa. 

its.  Tha  axtandad  sourcaa  of  raflaction  aro  dapandant  on  tho  p,B 
eomprossion  profilaa  in  such  a  way  that  tha  croas  darivativo  tarn 
can  combino  with  tha  right  hand  sida  to  prodnco  spatially  similar 
sonrca  tarms  whan  tha  Eg  fiald  is  primarily  inductiva  and  tha  com- 
prossion  is  spatially  homoganaous.  This  fact  is  oxploitad  indiraetly 
in  tha  naxt  saction. 

III.  Special  Cases 


Qna  can  restrict  tha  problem  to  quasi>static  fields  (drE  <  J)  and 
seek  spacializad  solutions  to  tha  momentum  conservation  equation  . 

In  such  cases  tha  Eg  fiald  is  a  proscribed  spaco-timo  function  and 
dominantly  inductiva  in  character,  axtarnal  electrical  rasarvoirs  sup¬ 
ply  any  naadod  currants,  and  tha  solution  defined  inside  some  radius 
can  ba  viewed  as  a  circuit  elamant. 

Using  the  conserved  flux  V<('’o)  =  dr  Sj(r)  to  formulate  the 
J  X  B  stress  term,  tha  momentum  aquation  can  ba  written 


DtUr  =  R 


ro  dR^  dR'  Rp 


where  Hj_{ro)  =  Be{ro)W ^{ro)  ,  and  c\{ro)  =  [dro^<\^ /^^Po 


If  oa«  has  tha  fnrthar  raatrietioaa  that  H±  ,  ara  indapandant 
of  fo  ,  aad  that  tha  iaitial  Beirg,)  is  iabaddad  such  that 
r~^  Cvt(ao)  .  a  coastaat,  thaa  a  apaeific  radial  daaaity 

profila  caa  liak  tha  tvo  taaaor  praaanra  eoapoaaata  —  tha  iaothar- 
aal  parpaadicalar  aaargy  aad  tha  rigid  rotor  eirenlatioa  oaorgy  —  vil 
tha  aaifora  Alfvaa  eonpraaaioa  profila.  Tha  iatrodnctioa  of  tha  tarn- 
porally  acalad  profila 


n(r,/) 


exp  — 


[  ( 


_r _ ^ 

<*0  A(f) 


r  +  ( 


A(0lf| 

2H± 


prodneas  a  liagla  aqaatioa  of  aotioa  for  A(f)  at  aay  radiai  <  r\ 
tha  radios  of  vaaishiag  carraat  daaaity.  Tha  total  carroat,  within 
r/(/)  =  A(/)  rj  .  for  N  ioas  par  ea.  is 


[k  .  i,5  ^  |i, 


whara  leennett  =  {2NmiH^yf^  ,  cs  =  (2  .  e  =  cs/c^  , 

5  =  ,  and  bg  =  •  For  a  given  current,  the 

effective  electron  temperature  la  Tx  =:  3.47 •  ^ 

Tj_  =  0.994  [kev]  at  Iiotai  =  1*26  MA  ,  N  =  5.0- 10^®  cm~^  ,  e  =  0.1  . 


Such  a  typical  profila  sat  is  shown  in  Fig.  2. 
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Fig.  2  Typical  Solution  Profiles 

3.55- cm-’  O  .Be  :  1.45  •  10<  Gs  X 
4.2stV/cm  □  .Ur  :  9.36  •  10®  cm /sec 


formalizing  tiio  timo  by  tbs  soand  spood  associatod  vitb  H±  .  lot 
ts  =  Oo/cs  ,  tsdr  =  dt  .  Tho  scale  factor  aqaation  of  motion  in  "sound 
transit  timo  units"  is  than 


A 


—  f— 

AM  Um 


with  initial  conditions  A^  =  .  ^(^o)  =  1  ■  This  aquation 

has  an  snargy  integral, 


h  =  I  +  Aj/2 


( 


1  + 


A  InX^ 


For  £  <  1,  Ao  =  0  there  is  always  an  interior  turning  point  and 

oscillation  inside  A  =  1  ;  for  £  >  1,  Ao  =  0  there  is  an  exte¬ 

rior  turning  point  and  oscillation  outside  A  =  1  .  Small  values  of 
£  imply  implosion  at  velocities  higher  than  C5  and  deeper  turning 
points  arrived  at  more  quickly.  Soma  trajectories  and  turning  points 
are  shown  in  Fig.  3,  labeled  by  various  Ao  ,  e  combinations. 


Such  a  gyrokinatic  oscillator  represents  a  special  case  of  the 
propagation  problem  for  which: 


^m«n(  ^01^  ) 

A,  =  (-0.0,  -5.0] 
c  =  [0.05, 0.95] 


a.  th*  Alfvan  spaed  is  linear  in  space  and  invariant  on  a  Lagrangian 
trajectory,  and 

b.  tlia  larger  total  current  cases  are  those  for  which  the  electron 
population  is  hotter  and  the  Alfven  speed  is  also  large  compare  to 
the  sound  spaed,  and 

c.  the  electromagnetic  energy  propagation  is  slowed  further,  below 
the  Alfvsn  speed,  by  the  tensor  pressure  of  the  particle  population 
and  suffers  a  complete  reflection  in  the  expansion  phase  of  the  mo¬ 
tion. 

The  present  solution  contrasts  with  a  comparable  MHD  solution  due 
to  Felber^  in  that  the  implosion  remains  hollow  on  axis,  recoiling 
from  the  "centrifugal  potential"  at  soma  nonzero  radius.  In  addition, 
the  tensor  pressure  requires  no  "completing"  equation  of  state  assump¬ 
tion,  being  specified  through  observables  at  soma  initial  time.  The 
currant  density  is  also  constituted  in  a  fully  self-consistent  man¬ 
ner,  requiring  no  "infinite  conductivity"  assumption  in  a  linear  Ohm's 
law. 

The  energy  absorbed  over  one  quarter  cycle  .^(e,  Ag)  ,  i.a.  down 
to  the  turning  point,  can  be  conveniently  written  in  terms  of  the  to¬ 
tal  current  at  any  typical  value  €o  and  the  logarithm  of  the  turn¬ 
ing  point  value  A^tn 


and  H{e,Xo)  {f-)^ 


is  shown  below  for  hotai  =  1-26  MA,  m,  =  12  amu, 


/V.  Implications  for  DNA  Objectives 


If  gyrokinatic  cnrrant  densities  approaching  those  of  the  salf- 
siailar  oscillator  are  actually  present  in  PES  diodes,  than  one  mnst 
suspect  a  significant  loss  of  driver  current  to  the  main  pinch.  The 
power  absorbed  in  the  corona  plasma  may  be  mnch  more  modest ,  since  it 
depends  on  the  difference  of  incoming  and  outgoing  Eg  field  com¬ 
ponents,  but  the  possibility  of  significant  reflected  power  from  the 
corona  cannot  be  discounted.  However  this  is  the  worst  case  scenario. 

It  is  more  likely  that  gyrokinatic  populations  are  difficult  to  ob¬ 
tain  in  the  clean  form  required  by  our  simple  solution.  They  may  tend 
to  occur  only  once  in  a  while  and  participate  in  the  many  factors  which 
produce  the  shot  to  shot  variation  in  PES  diodes.  It  is  clear  from 
the  theory  to  date  that  a  serious  experimental  effort  is  warranted  to 
fix  the  early  peripheral  conditions  which  help  define  the  gyrokinatic 
corona.  Experiments  of  particular  importance  would  be  early  time  neu¬ 
tral  gas  density  measurements  as  (r,  2)  profiles  and  spacetime  resolved 
magnetic  field  maasuramants  emphasizing  the  period  of  run-down  and 
stagnation.  Since  the  radial  extant  of  the  gyrokinatic  corona  can  be 
on  the  order  of  centimeters,  spectroscopic  (Zeeman  effect)  measure¬ 
ments  based  on  seeded  trace  impurity  concentrations  may  prove  viable. 
Such  work  can  be  guided  by  more  general  theoretical  calculations  of 
the  power  flow  which  are  made  possible  by  the  progress  achieved  so  far. 

Even  mors  intriguing  is  the  issue  of  microstability  for  the  Vlasov 
solutions  examined  here.  Having  fixed  only  the  lowest  velocity  mo¬ 
ments  in  definins  the  current  density,  a  infinite  sequence  of  higher 
momenta  remains  available  and  the  quasilinear  spectrum  of  unstable 
modes  which  might  be  excited  remains  to  be  determined  through  these 
observables.  Insofar  as  the  gyrokinetic  populations  have  been  shown 


capable  of  suBtaiaing  high  currents,  these  populations  may  veil  prove 
to  be  useful  as  sivitches  if  the  conditioning  of  the  higher  moments  can 
produce  unstable  plasmas  which  vill  interrupt  those  currents  in  a  pre¬ 
dictable  time  frame.  The  present  formulation  is  perhaps  the  first  the¬ 
oretical  development  suited  to  the  study  of  such  phenomena  in  a  sim¬ 
ple,  yet  geometrically  relevant  context. 
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F .  Eroding  Foil  Switch  Model  for  Mass  Flow  Switches 

I.  THE  EPS  MODEL  AHD  FIELD  PENETRATION 

Ilia  Eroding  Foil  Switch  (EFS)  oodal  waa  originally  foranlatad^  as 
a  zaro  diaansional  pietnra  to  fit  nnaxplainad  axpariaantal  raanlts. 

It  was  darivad,  not  from  first  prineiplas  bnt  rathar  ampirically,  from 
study  of  axparimantal  rasults  and  physically  raasonabla  assumptions. 
Prsliminary  rasults  of  both  thaory  and  numarical  simulation  for  tha 
magnatic  diffusion  procass  in  thin  foils  hava  shown  tha  basic  concepts 
of  tha  modal  to  ba  accurate  and  a  more  thorough  investigation  is  on¬ 
going. 

Vhilo  tha  accelaration  of  foils  in  a  coaxial  gsometry  is  truly  a  two 
dimensional  problem  due  to  variation  of  tha  magnetic  field  strength 
with  radius,  it  is  claar  from  axarcisaa  with  the  simulation  coda 
iXFQIL^  that  proper  "mass  grading*  can  produce  a  stably  acealaratad, 
essentially  one  dimansional  motion.  This  suggests  that  models  for  tha 
kinematics  and  fiald  dynamics  which  ara  restricted  to  the  axial  do¬ 
main  can  remain  viabla  analysis  tools. 

In  events  like  tha  "foil  switch"’  and  "magnetic  gata"^  axparimants 
at  LAEL  tha  foil  behavior  indicates  that  soma  substantial  mass  loss 
is  a  prims  feature.  For  both  axparimants  the  object  is  to  accalar- 
ata  a  foil,  using  Lorantz  forces,  past  a  gap  in  the  transmission  line 
where  switching  can  occur.  The  foil  switch  sxparimant  is  dona  in  a 
coaxial  transmission  line  with  one  switch  gap,  while  the  magnatic  gate 
sxparimant  occurs  in  a  linear  striplina  with  three  switching  gvps.  In 
Fig.l,  2,  and  3  tha  gaomatry,  currant  and  foil  position  ara  summarized 
to  indicate  that  as  little  as  five  percent  of  tha  original  mass  may 
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partieipat*  ia  tha  final  eoamatation  of  tha  carrant.  Tha  awiteh  sig¬ 
nal  itaalf  was  a  liaaar  risa  ovar  10~^  sac,  oaa  taath  tha  charaetar- 
istie  tiaa  asaoeiatad  with  tha  foil  aquation  of  notion  with  tha  mass 
fixad  at  its  original  valna.  This  indicatas  a  broad  cnrrant  shaath, 
parhapa  a  cantinatar  wida,  as  tha  nachanisa  of  currant  transfar  to  tha 
"load*  cavity.  Tha  loat  naas  was  also  hot,  optical  pyromatry  indi- 
catad  a  foil  tamparatura  of  4000  ^C.  Tha  aass  loss  was  smooth  in  tiaa 
rathar  than  abrupt  and  a  fair  fraction  of  tha  foil  had  apparantly  va- 
poriaad. 

Sinca  tha  foil  aass  had  baan  choaan  to  insnra  much  lowar  taapara- 
tnras  undar  tha  assumption  that  tha  cnrrant  would  ba  uniformly  dis- 
tributad,  it  saaaa  that  soma  localization  of  tha  currant  is  raspon- 
sibla  for  both  tha  high  taaparatnraa  and  tha  rapid,  nonnniform  accal- 
aration  of  tha  switch  mass.  Concantration  of  tha  currant  in  a  thin 
snrfaca  layar  would  prodnca  aarly  haating  and  mass  arosion.  Tha  highar 
tamparatura  incraasaa  tha  local  rasistanca  and  forcas  tha  currant  dan- 
sity  to  propagata  daapar  into  tha  natal.  As  tha  procass  continual  tha 
haatad  portion  of  tha  switch  aithar  malts  or  vaporizas,  but  tha  point 
of  maximum  axial  straas  is  always  in  front  of  tha  high  tamparatura 
zona.  Soma  of  tha  liquid  phasa  and  all  of  tha  vapor  phasa  is  than  laft 
bahind  and  tha  axial  kinamatics  is  datarainad  by  tha  aass  ahaad  of  tha 
Jr  X  Be  strass  front . 


A.  Eroding  Slab  Picture 

Eagardlass  of  tha  datails  in  tha  mass  loss  process ,  tha  rata  of  mass 
loss  should  ba  proportional  to  tha  ohmic  haating  rata  and  thus  to  tha 
square  of  tha  magnetic  field  on  tha  upstraam  boundary  of  tha  foil.  As 
a  rough  cut  at  the  erosion  procass,  one  may  than  let 


dtn 


1 
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where  )i  is  the  foil  mass  per  unit  area,  and  K  is  measured  in  [^/cm  sec 
baij  when  the  magnetic  field  pressure  is  measured  in  bars,  1  bar  =  1.0-10°  Gs 
Since  the  mass  graded  foil  equation  of  motion  for  V^(t)  can  be  integrated 
from  rest  to  yield 

„V,  =  W[l)  =  j‘dt,  , 

tha  axial  valoeity  of  tha  arodiag  foil  can  ba  axpraaaad  in  tarma  of 
tba  valoeity  Ve(^)  =  W{t)/iio  of  a  fixad  mass  foil.  viz. 


'  (1  -  BT  V.(<))  ’ 

zltb  fio  tba  initial  foil  maas.  In  tbaaa  axpraaaiona  tba  fiald 
valna  SBg  eonnotaa  a  fiald  praaanra  difference  aeroaa  tba  alab.  Aa 

maaa 

arodaa  tbia  praaanra  differanea  ia  appliad  daapar  in  tba  alab;  tbia 
rough  model  aaaumaa  that  arodad  maaa  ia  highly  raaiativa  —  tba  orig¬ 
inal  aurfaea  SB^  ia  Juat  ahiftad  to  the  interior.  In  order  for  tha 
modal  to  retain  aalf-conaiatancy,  tha  factor  K  muat  ba  indapandant 
of  radiua.  In  any  firat  principlaa  davalopment  of  thia  theory  thara- 
fora  one  muat  head  whatever  eonatrainta  are  needed  to  maice  thia  hap¬ 
pen. 

Tha  factor  K  haa  bean  evaluated  in  three  vaya  for  tha  paramatara 
of  tha  LANL  experiment.  Tha  firat  ia  through  a  finite  difference  cal¬ 
culation  of  tha  aroaion  proceaa;  tha  aacond,  through  a  aalf-aimilar 
aolntion  for  tha  fiald  penetration;  the  third,  by  examination  of  the 
experimental  kinamatica.  All  three  evaluationa  are  in  rough  agree- 


B.  Magaatie  Diffnaiea  aad  Salf-Siailar  Daaeriptions 

Tlia  aagaatic  dillasioa  aqaatioa  la  aatarally  poaad  la  taraa  of  di- 
aaaaioalaaa  variablaa  r  =  tjto  ,  i  =  z//o  «  •  Tha  axial 

valocity  of  tlia  foil  ia  thaa  =  Vzlc  .  aad  tha  affaetiva  raaia- 
ti7it7  ia  i\  =  I  f  toff  .  Tbia  raaiativity  ia  gaaarally  a  faactioa 
of  tba  iataraal  aaargy  par  aalt  aaaa  of  tba  awitcb  aatarial.  Ia  tba 
followiag  davalopmaat  tha  baaic  aiaalatioa  nait  ia  a  riag  at  aoaa  ia- 
tarmadiata  radina  to  ia  a  coaxial  wavagaida,  with  a  aaaa  hiring  • 
Tha  alaetromagaatie  powar  dapoaitioa  ia  thia  riag  ia  tha  apatial  ia- 
tagral  of  radial  carraat  daaaity  aqaarad  aad  tha  local  raaiativity, 
viz. 


DrQ  =  »(<3(?,  r)/M)  S.l’  . 

whara  hg  =  ly+j  -  ly  tha  axial  call  aiza,  Ao  =  ('‘f  +  A?)  “  tha 
croaa  aactioaal  araa  factor,  aad  DrQ  ia  tha  aatarial  darivativa  with 
raapact  to  tha  acala  tiaa  of  tha  iataraal  aaargy  par  ring.  latagra* 
tion  of  thia  aqaatioa,  ovar  aay  tiaa  doaain  oa  which  tha  axial  gra- 
diaat  of  Bg  ia  providad,  allowa  tha  aalf-coaaiataat  evalnatioa  of 
tha  raaiativity  aa  ri{Q{z,T)lhirin9)  ■  Variona  raaiativity  aodela  ara 
availabla  and  the  praaant  iapleaantation  ariaaa  froa  a  fit  to  LA2fL  and 
AFWL  data  for  Alaainna^  of  the  fora  7  =s  [l+.4(Q/.'V/)+B((?/iV/')^]  . 

Tnrning  to  tha  aa,^atic  field,  ona  nay  write  tha  diffnaion  aquation 
with  tha  aaaa  aatarial  darivativa  operator  applied  to  Bg 

DrBg  =  ^  di  \  n{Q[z,T)lM)  d-zBg]  -  Bg  dg^g. 

47r 
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Ilia  eoavaetiva  tara  iavolTing  tlia  axial  gradiaat  of  can  ba  ig¬ 

nored  in  tha  liaplaat  EFS  application  baeansa  tba  material  ia  aaanmad 
ineompraaaibla  until  it  ia  ia  tba  vapor  phaao  and  long  loot  to  tba  maaa 
being  coavaetad. 

Ibia  formulation -will  laava  tba  eolation  invariant  vitb  tba  radiaa 

» 

To  provided  tba  ring  maaa  ia  graded  aneb  that  tba  off active  maaa  den- 
aity  acalea  a  {Rgeaul^a)^  bba  driving  boundary  valna  ia  acaled 

[Rie^iel^o)  •  Tba  axial  acalaa  muat  ba  praaarvad  in  tba  tranafor- 
mation,  bovavar,  ao  tba  EFS  theory  ia  limited,  rigoroaaly,  to  tboaa 
maaa  graded  caaaa  wbara  tba  avitcb  maaa  denaity,  ratbar  than  tba  ax¬ 
ial  foil  tbicknaaa,  la  conditioned. 

With  tba  above  caveat  in  mind,  to  aaaaaa  tba  EFS  maaa  aroaioa  coa- 
ataat  K  oaa  naad  only  anpply  an  apatraam  boundary  condition  for  tba 
fiald  Be  ,  togatbar  vitb  a  avitcb  maaa  and  radiua  vbicb  are  conaia- 
tant  vitb  tba  experiment.  Tba  raaulta  are  illnatratad  in  Fig.  4. 

rolloving  tba  currant  daanity  peak  acroaa  figuraa  aucb  aa  tbaaa  oaa 
can  aaaaaa  tba  maaa  uncovered  by  ita  prograaa  and  tbna  aatabliab  a  fac¬ 
tor  K,  at  tba  foil  maaa  and  for  tba  driver  current  wavaform  cboaaa 
in  any  modal  aolution.  At  praaant  tbia  analyaia  indicataa  a  value 
Kpde  =  27  [gm/cm^  aac  bar] ,  but  planned  modal  refiaamanta  and  mora 
detailed  exploration  of  tba  initial  eonditiona  will  probably  reduce 
tbia  raault. 

An  altarnata  path  to  aucb  an  aatimata  liaa  in  aalf-aimilar  or  pro- 
groaaive  vava  aolutiona  to  tba  diffuaion  and  beating  aquationa.  Hera 
of  conraa  tba  apacializad  boundary  eonditiona  will  not  accomodate  gen¬ 
eral  driver  waveforma,  but  tba  aaaantial  alamanta  of  tba  penetration 
pbyaica  can  ba  diaplayed.  Let  A  =  tioAjM  v  4.8-10~^^  [aac  bar], 
uaing  tba  roaiativity  axpanaion  above.  In  tba  caaa  ,3g=0,  conaidar 
tba  diffuaion  aquation  in  tba  region  wbara  dgJr  =  d^Be  =  0  .  At 
tbia  point  tba  diffuaion  aquation  ia  a  aimpla  advactiva  relation,  viz. 


—  dgti  dzBe 


!  drBe 

\ 

t 


and  any  prograaaiva  wava  lolntioa  will  propagata  at  apaad  Vb  —  dzH- 
Uaiag  tha  haatiag  aquation  with  this  simpla  fiald  aquation  providas 
a  local  spatial  constraint  on 


dm 


A  n 

4r 


1/3 


dzBg. 


It  this  aquation  is  intagratad  ovar  an  intarval  6z  back  into  tha 
aagnatic  fiald  wava,  ona  can  ralata  ti  to  tha  changa  in  aagnatic  fiald 

A  S 

prassura,  obtaining  q  =  *■  .  Wharaupon  a  first  ordar  asti- 

aata  of  tha  gradiant  in  q  and  soma  simplification  providas  tha  apaad 


Vb 


j_A_ 

^32  6z 


Evaluating  in  cgs  units  and  racalling  that  dtfi  -  --K  SBjjSir  s 
—PAi  I  hha  "salf-similar"  arosion  factor  bacomas 


K,z  =  =  21  [gntjcm^  sec  6ar]. 

43r  0  z 

II.  The  EFS  Model  and  Foil  Kinematics 

As  tha  mass  erosion  proceeds  through  the  slab,  a  final  state  obtains 
wharsin  tha  foil  must  vaporize  to  a  slightly  ionized,  snowplowed 
plasma.  At  this  point  the  remaining  mass  will  convact  at  tha  local 
Alfvan 

shock  speed®  .  Using  tha  value  =  3.0-10®  cm/sac 

from  tha  LANL  switch  axparimantal  paramatars  and  assuming  a  changa  from 
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•roaion  to  eoavoetioa  at  a  tiaa  ,  tha  aqaatioa  to  fit  for  tha  ax' 
pariaaatal  kiaaaatics  caa  ba  writtaa 


(i-K  d'fl, )/»<.)  ■ 


By  aatcbiag  valoeitiaa  at  ts  from  tba  tvo  modal  phaaaa,  oaa  may  as- 
timata  ts  aad  tbaa  aafold  K  throagb  tba  match  coaditioa  itsalf, 


Mo  _  ir  - 
W{tE)  Be  ' 


For  tba  LiSL  svitcb  parametars  oaa  fiads  tE  =  3.9  ^aae  aad 
K  =  IQ  [gm/cm^  sac  bar] ,  somawhat  balow  tba  two  tbaoratieal  as- 
timataa . 

III.  Implications  for  DNA  Objectives 

It  tba  prasaat  laval  of  sopbisticatioa,  tba  EFS  modal  appaars  to  ax- 
plaia  maay  basic  faataras  of  mass  flow  switcbas  which  origiaata  tba 
switcbiag  actioa  with  a  solid  phasa  coadaetor.  If  a  wall  dafiaad  ax- 
parlmaat  is  fisldad,  tba  datailad  applicatioa  aad  varificatioa  of  tba 
modal  will,  if  saccaasfal,  provida  a  vary  asafal  * cal i brat ioa”  of  tba 
mass  orosioa  procsss.  A  asafal  axparimsat  would  iavolva  tbraa  crit¬ 
ical  alamsats: 

I.  Axial  rssolatioa  of  tba  locatioa  of  tba  curraat  cbaaaal,  aad 
perhaps  the  width  as  wall,  as  a  faactioa  of  tiaa.  Tbs  obvious  method 
would  ba  saquaatial  curraat  moaitors  aloag  either  coaxial  coaduc- 
tor,  but  tba  short  laagtb  of  travel  may  prasaat  something  of  a  cbal- 

laaga . 

It.  Temporal  resolutioa  of  the  voltage,  measured  downstream  across 
tba  switch,  during  tba  conduction  (run-down)  phasa.  This  acta  as 
a  diagnostic  of  motional  impedance  and  field  penetration  through  the 
foil.  Naturally,  it  suggests  tba  use  of  an  open  circuit  termina¬ 
tion  on  tba  coaxial  line  in  which  tba  foil  travels. 


ttt.  Temporal  roaolatioa  of  tho  apatraam  foil  tomporatara  and  tlia  fell 
blovoff  material,  with  particular  attaation  to  an  aatimata  of  the 
siaa  of  any  liquid  droplata  praaant  in  tha  eroded  material. 

Once  thoroughly  tasted,  tha  EFS  modal  can  be  incorporated  in  more  gen¬ 
eral  switch  models,  such  as  AIFQIL,  and  will  thus  improve  the  physics 
content  of  such  codas. 
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Figure  1.  A  Diagram  of  the  LANL  Foil  Switch  Experiment. 
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Figure  2.  Plot  of  normalized  switched  current  (switched 
current  divided  by  bank  current)  for  Los  Alamos  foil  switch 
experiment. 
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Figure  3.  A  plot  of  foil  position  versus  time  from  1) 
experiment  data  (squares)  2)  F  s  ma  calculations  using  the 
entire  foil  mass  (dotted  curve).  3)  F  s  ma  calculations  using 
5%  of  the  original  foil  mass  (solid  curve). 
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Fig.  4.  EFS  field  penetration  behavior. 
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G.  Proton  Seam  Deposition,  Heating,  and  Radiation  from  an  Fe  Plasma 


I.  INTRODUCTION 


Recent  advances  in  pulsed  power  and  diode  technologies  nave  made  it 
possible  to  deliver  intense  proton  beams  onto  a  target  with  high 
efficiency.  Thus,  the  interaction  of  energetic  charged  particle  beams  with 
different  materials  has  become  an  active  field  of  research.  A  significant 
application  for  this  research  is  the  use  of  ion  beams  for  x-ray 
generation.  Other  applications  include  and  weapons  lethality  and 
vulnerability  studies. 

The  interaction  of  a  monoenergetic  one  MeV  proton  beam  with  a  planar 

1  --j 

aluminum  slab  of  15uni  thickness  has  been  previously  investigated.  Tnis 
beam  was  assumed  to  consist  of  a  square-shaped  pulse  of  10  nanosecond 
duration  with  a  flux  of  10“^  protons  cm  sec  .  The  beam  intensity  on 
target  is  1.6X10^^  W/cm^,  which  is  comparable  to  intensities  available  with 
current  pulsed-power  generators. 

In  the  present  study,  the  interaction  of  this  same  proton  beam  with  a 

planar  8.0-um-thick  iron  slab  is  treated.  This  distance  approximately 

u 

corresponds  to  the  range  of  a  one  MeV  proton  in  cold  solid  density  Fe. 

The  Fe  atomic  model  used  in  this  calculation  did  not  have  excited 
level  structure  in  ionisation  stages  below  Fe  XVI.  Only  ground  states  were 
carried  for  stages  I  through  XV.  While  this  lack  of  level  structure  and 
the  resulting  omission  of  lines  from  the  radiation  transport  scheme  can 
modify  the  energy  balance  in  regimes  where  these  lower  stages  are  highly 
populated,  it  is  estimated  that  the  effect  is  less  than  a  factor  of  two  on 
the  energy  lost  due  to  radiation. 

The  scheme  for  treating  the  proton  beam-iron  slao  interaction  is 
fundamentally  the  same  as  that  described  in  Refs.  '-3  for  treating  the 
proton  beam-Al  slab  interaction.  The  hydrodynamic  response  is  treated  via 
a  one-dimensional  model  with  a  sliding-zone  version  of  flux  corrected 
transport.  A  special  gridding  algorithm  moves  zones  in  a  Lagrangian 
fashion  and  adjusts  the  mesh  in  order  to  resolve  steep  gradients  in  the 
flow.  A  single  temperature  is  assumed.  A  collisional  radiative 

equilibrium  treatment  of  the  ionization  balance  is  used.  Radiation 


transport  is  treated  with  a  hybrid  scheme  which  uses  multi  frequency  ray 
tracing  for  the  continua,  and  a  probability-of-escape  method  for  lines. 
The  radiation  transport  is  coupled  self -consistently  into  the  ionization 
dynamics.  The  proton  stopping  power  includes  contributions  from  bound 
atomic  electrons,  free  plasma  electrons,  and  plasma  ions. 

A  major  difference  between  this  work  with  Fe  as  the  target  material 
and  the  previous  work  with  Al  as  the  targ*et  material  is  in  the  treatment  of 
the  bound  electron  stopping  power.  For  the  proton-Al"^^  interaction,  where 
q  is  the  ionic  charge,  detailed  stopping  power  cross  sections  are 
available;^  for  the  proton-Fe*'^  interaction,  no  such  data  are  currently 

C 

available.  In  a  previous  simulation  of  a  proton  beam-Fe  slab  target 

T  8 

interaction,  a  modified  local-oscillator-model  (LOM)'’  was  used  to 
estimate  stopping  power  due  to  bound  electrons.  Figure  1  (taken  from  Ref. 
3)  shows  an  LOM  calculation  of  proton  stopping  power  in  neutral,  solid 

M 

density  Fe  in  comparison  with  the  Andersen-Ziegler  data  for  this  case.  In 
this  case,  the  LOM  gives  good  results  at  low  and  high  energies,  but  it 
underestimates  the  peak  and  shifts  it  to  lower  energy.  This  fact,  coupled 
with  some  problems  with  the  line  spectra  (which  have  since  been  resolved) 
in  the  previous  treatment,^  led  to  a  second  attempt  to  treat  the  proton 
beam-Fe  target  interaction. 

In  the  present  treatment,  an  analytic  fit  to  the  Andersen-Ziegler 
proton-neutral  Fe  stopping  power  was  obtained.  Then,  stopping  power  for  Fe 
ions  was  scaled  from  the  neutral  Fe  stopping  power  by  the  number  of  atomic 
electrons  on  the  ion;  i.e., 

S^(p  ^  Fe^'^)  -  ((  Z-Z  )  /  Z)  »S^(p  *  Fe*°)  , 

where  is  the  bound  electron  stopping  cross  section,  Z  is  the  atomic 
number,  and  Z  is  the  average  charge  of  the  target  Fe  ions. 

There  are  some  difficulties  with  this  stopping  power  model.  One  is 
that  the  peak  stopping  power  occurs  at  the  same  proton  energy  for  all 
ions.  In  the  detailed  Al*^  stopping  power  calculations^,  the  peak  shifted 
to  higher  energy  as  q  increased.  Another  difficulty  is  that  Sq.  ('  ) 
probably  overestimates  stopping  power  due  to  the  ions.  The  detailed  Al*^ 
calculations  show  the  ratio  of  ion  to  neutral  stopping  power  to  be  lower 
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than  that  indicated  by  Eq.  (i  ).  However,  due  to  the  lack  of  more  detailed 
information,  some  assumptions  regarding  stopping  power  must  be  made  to 
treat  the  problem. 

Some  results  of  the  most  recent  proton  beam-Fe  slab  target  interaction 
simulation  will  now  be  given  and  discussed.  Figure  2  shows  the  energy 
history  and  partitioning  in  the  beam  target  interaction.  The  effect  of  the 
beam  cutoff  at  10*  nsec  can  be  clearly  seen  as  the  total  plasma  energy 
decreases  after  this  time.  This  is  due  to  thermal  (internal)  energy  losses 
as  the  plasma  radiates  and  to  potential  (ionization)  energy  as  the  ions 
recombine.  The  radiated  energy  is  negligible  at  early  times,  but  becomes  a 
significant  factor  at  late  times;  at  13.2  nsec,  it  represents  about  45  per 
cent  of  the  deposited  energy. 

Figure  3  gives  temperature  and  density  profiles  at  different  times  in 
the  plasma  evolution.  Generally  the  plasma  gets  wider  and  less  dense  as 
time  goes  on.  However,  the  temperature  appears  to  be  relatively  constant 
over  most  of  the  plasma  at  later  times.  This  is  reflected  in  the  fact  that 
the  internal  energy  curve  is  fairly  constant  during  the  later  times  until 
the  beam  is  turned  off.  The  13.248  nsec  temperature  profile  is  slightly 
cooler  than  the  10  nsec  profile. 

Figure  4a  displays  the  temperature  and  density  profile  at  1.59  nsec 
from  the  previous  proton  beanrFe  target  treatment.^  Although  the  plasma 
widths  are  approximately  the  same,  the  rear  surface  temperature  is  lower 
and  the  density  higher  in  Fig.  3(a)  at  1  .569  nsec.  The  front  surface 
temperatures  are  about  the  saaw  in  the  two  cases,  but  the  present  treatment 
yields  a  much  lower  front  surface  density.  The  differences  in  these 
profiles  is  due  to  the  fact  that  Eq.  (1)  gives  larger  stopping  powers  than 
the  L.0M;  hence  the  beam  does  not  penetrate  as  far  into  the  target  and 
deposits  its  energy  nearer  the  front. 

Figure  4b  shows  temperature  and  density  profiles  at  9.75  nsec  in  the 
previous  treatment  of  an  Fe  target.^  Comparison  with  the  10  nsec  profile  in 
Fig.  3  shows  that  the  previous  treatment  produced  much  more  uniform  plasma 
profiles  near  beam  shutoff,  while  the  present  treatment  produces  distinct 
variations  in  temperature  and  density  near  the  rear.  The  front  surface 
values  are  about  the  same  in  the  two  treatments.  These  differences  are 
again  due  to  the  greater  penetration  allowed  by  the  COM  stopping  power 
model. 
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Figure  5  shows  the  time  behavior  of  the  density,  temperature,  and 
average  charge  Z  at  the  front  and  rear  surfaces.  At  the  front  side, 
temperature  and  Z  rise  very  rapidly  and  then  maintain  fairly  steady  values 
afterward;  the  density  drops  rapidly  at  first,  then  drops  more  slowly  to 
about  0.01  gm/cm^  around  13  nsec.  At  the  rear  surface,  density  and 
temperature  change  very  rapidly  initially,  then  the  variation  proceeds  more 
slowly.  They  do  not  achieve  steady  values  however.  The  Z  curve  follows 
the  temperature  curve.  The  peaks  in  the  curve  at  6  nsec  are  due  to  a  shock 
arriving  at  the  rear. 

In  Fig.  6  are  displayed  the  calculated  front  and  rear  spectra  at  the 
same  times  as  the  plasma  profiles  in  Fig.  3.  The  lack  of  level  structure 
in  ion  stages  below  Fe  XVI  is  clearly  shown  by  the  scarcity  of  lines  at 
1  .569  nsec.  As  might  be  expected  from  the  density  and  temperature  plots, 
most  of  the  radiation  is  lost  from  the  front  of  the  plasma;  most  of  this 
radiation  comes  from  zones  near  the  front  of  the  plasma.  For  the  rear 
spectra,  the  least  energetic  radiation  arises  from  zones  near  the  rear. 
The  increased  emissions  near  0.7  keV  and  around  2  to  3  keV  consists  of 
contributions  from  Interior  zones  which  escape  through  opacity  "windows". 
At  10  nsec  and  13.2U8  nsec,  the  radiation  rises  due  to  contributions  from 
interior  zones.  These  bumps  in  the  rear  spectra  coincide  with  minima  in 
the  optical  depth  curves  allowing  radiation  emitted  in  the  interior  to 
"shine  through"  to  the  rear  aide.  To  a  lesser-  extent,  the  same  phenomenon 
accounts  for  the  relative  maxima  and  s'.iima  in  the  front  spectra.  In  the 
previous  Fe  target  simulation,®  the  rear  surface  emission  was  much 
higher.  In  that  case,  the  front  and  rear  spectra  at  9.76  nsec  and  12  nsec 
were  very  nearly  equal.  These  differences  are  presumably  due  to  the 
differences  in  stopping  power  predicted  by  the  LOM  and  used  in  this  model. 

The  results  of  this  study  and  its  comparison  with  a  previous 
treatment®  of  the  same  problem  reinforces  the  necessity  of  having  good 
stopping  power  data  in  ion  beam-target  interaction  studies,  particularly  in 
higher  Z  targets  where  bound  electron  stopping  power  is  significant  over  a 
wide  range  of  temperatures  and  densities.  In  low  Z  ions  where  free 
electron  stopping  becomes  dorainent  at  low  Z,  accurate  bound  electron 
stopping  power  becomes  less  critical  at  later  times. 
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The  Fe  model  used  here  has  no  structure  In  ion  stages  below  Fe  XVI; 
hence,  radiation  transport  is  affected  at  all  stages  of  this  simulation, 
since  the  maximum  Z  attained  was  about  16.  Radiation  losses  are  hence 
underestimated,  particularly  near  the  rear  surface,  but  it  is  estimated 
that  the  error  is  less  than  a  factor  of  two.  It  was  found  that  radiation 
transport  and  losses  are  a  significant  factor  in  beam-target  interactions. 

This  is  a  one  dimensional  treatment;  hence,  field  effects  and  energy 
flow  perpendicular  to  the  beam  direction  are  omitted.  These  factors  could 
change  the  profile  shapes  and  magnitudes,  and  the  projectile  ranges  in  the 
plasma,  and  hence  alter  the  results  presented  here. 


References 

1.  J.  E.  Rogerson,  R.  W.  Clark,  and  J.  Davis,  Phys.  Rev.  A  21t3323  (1985). 

2.  J.  E.  Rogerson,  R.  W.  Clark,  and  J.  Davis,  MRL  Memorandum  Report  5^82 
(1984). 

3.  Plasma  Radiation  Branch,  PPD,  NRL  Memorandum  Report  5590  (1985). 

4.  H.  H.  Andersen  and  J.  F.  Ziegler,  Hydrogen-Stopping  Powers  and  Ranges 
in  all  Elements,  Pergamon  Press  C977). 

5.  S.  J.  McGuire,  J.  M.  Peek,  and  L.  C.  Pltchford,  Phys.  Rev.  A  26, 13' 3 

(T982) . 

0.  Plasma  Radiation  Branch,  PPD,  MRL,  DMA  Progress  Report  (5  July  1985). 

7.  R.  X.  Mesbet  and  J.  F,  Ziegler,  Appl,  Phys.  Lett.  21.310  (1977). 

I 

3.  J.  S.  Rogerson,  NRL  Memorandum  Report  4485  (198")  ADA  097772. 

I 

I 


92 


H.  Won-LTE  Electron  Dlatributlona  In  Electrical  Jischarges 

I.  INTRODUCTION 

The  electron  energy  dietrlbution  function  is  the  connecting  link 
between  the  electrical  and  radiative  behavior  of  electrical  discharges: 
amely,  it  determines  both  the  strength  of  the  discharge's  electrical  and 
thermal  conductivities  as  well  as  the  strength  of  the  ionization  and 
excitation  rates,  and  therefore  the  radiation  loss  rates.  Thus,  it  is 
essential  that  electron  distribution  functions  be  calculated  for  pulse 
power  x-ray  machines,  which  should  be  designed  to  be  as  energy  efficient  as 
is  physically  possible  in  converting  the  heated  free  electron  energy 
distributions  to  kilovolt  x-rays.  Without  a  good  knowledge  of  the  dynamics 
of  the  electron  distribution  function,  the  full  potentialities  of  pulse 
power  x-ray  generators  cannot  be  well  assessed. 

From  the  electrical  and  thermal  conductivities,  one  calculates  the 
current  distributions  (and  thus  the  details  of  the  Implosion  dynamics). 
Moreover,  it  has  long  been  appreciated  that  the  conduction  properties  of 
pulse  power  loads  are  anomalous,  i.e.  they  are  not  understood  in  terms  of 
classically  calculated  Spitzer  conductivities.  Within  the  DOE  ICF 

community,  flux  correctors  have  long  been  in  use  to  sharply  reduce 
classical  heat  conduction  in  steep  temperature  and  density  gradients. 
Experience  has  also  taught  this  community  that  the  values  of  the  flux 
corrector  can  change  markedly  depending  on  experimental  conditions.  Its 
use  also  contributes  to  the  poor  agreement  that  is  usually  obtained  in 
matching  theory  with  detailed  experimental  x-ray  data. 

II .  PAST  TREATMENTS 

There  are  several  inconsistent  features  to  past  treatments  of  the 

electron  distribution  function  problem.  For  ionization  calculations,  it  is 

generally  assumed  that  the  electrons  have  Maxwellian  distributions  f‘(E), 

which  is  equivalent  to  assuming  that  the  free  electrons  are  in  local 

2 

thermodynamic  equilibrium  (LTE).  E  ■  1/2  m  v  is  the  electron  kinetic 
energy.  In  order  to  compute  electrical  and  thermal  conduct! vites,  however, 
one  must  assume  from  the  beginning  that  the  electron  distribution  function 
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f  is  non-Maxwellian  and  that  it  depends  doth  on  the  electron  energy  E  and 
on  the  cosine  of  the  angle  u  that  the  electron  notion  has  with  respect  to 
the  net  effects  of  electric  fields  and  plasma  gradients.  In  lowest  order 
approximation 


f(S,u)  -  f  (E)  *  u  f,(E) 

O  I 


f'j  is  then  determined  as  a  local  function  of  f^  under  the  often  false 
assumption  that  elastic  collisions  are  always  dominant.  Spltzer 
conductivity  is  calculated  in  this  manner  and  is  widely  used«  with  flux 
correctors,  in  hydrodynamics  codes.  However,  electrical  and  thermal 
conductivities  calculated  in  this  way, 

1.  do  not  treat  the  effect  of  inelastic  collisions  on  f  and  therefore 
apply  only  to  fully  ionized,  hydrogen Ic  plasmas 

2.  assume  Maxwellian  f^ 

3.  are  not  valid  at  high  electric  field  strengths  or  in  steep 
temperature  or  density  gradients 

4.  do  not  treat  transient  time  behavior,  i.e.  any  effects  due  to  the 
rapid  rise  in  electric  field. 

Both  electrical  and  thermal  conductivities  and  excitation  and 
ionization  rates  depend  primarily  on  the  shape  of  the  high  energy  tail  of 
f.  Thus,  the  energetics  of  pulse  power  x-ray  simulators  has  two  key 
features.  Hot  only  must  the  electrons  be  bulk  heated  but  this  energy  must 
flow  continuously  into  the  high  energy  tail  of  the  electron  distribution  to 
replenish  the  energy  that  is  lost  when  electrons  collide  inelastically  with 
ions  that  radiate  away  the  excitation  energy.  Energy  that  is  fed  into  the 
high  energy  tall  of  f  is  also  associated,  however,  with  the  problem  of 
runaway  electrons  that  decouple  from  the  plasma  altogether  and  lose  their 
energy  by  colliding  with  the  anode.  However,  r’unaway  behavior  has  only  been 


analyzed  in  the  absence  of  inelastic  collisions,  whose  cross-sections  do 
not  fall  off  as  rapidly  with  energy  as  do  elastic  cross-sections.  Hence, 
the  runaway  electron  problem  has  not  been  properly  analyzed  to  date  though 
it  clearly  bears  on  a  potentially  major  energy  inefficiency  mechanism. 

To  briefly  summarize  the  present  situation,  plasma  conductivities  and 
the  details  of  power  flow  are  both  interrelated  and  time  evolving 
characteristics  of  imploding  z-pinches  that  are  impossible  to 
experimentally  measure  directly.  They  must,  therefore,  be  self- 
consistently  calculated.  However,  in  the  past,  the  few  people  who  have 
calculated  the  change  in  ionization  level  caused  by  non-Maxwe Ilian 
electrons  have  failed  to  correlate  this  effect  with  changes  in  electrical 
conductivity.  Also  their  calculations,  were  carried  out  only  for  hydrogen 
plasmas  and  showed  large  reductions  in  the  number  of  high  energy 
electrons.  In  a  z-plnch,  these  electrons  are  needed  in  abundance  to  excite 
the  K-shell  of  ions,  thereby  producing  the  kilovolt  radiation.  On  the 
other  hand,  people  who  have  calculated  electrical  conductivities,  find  that 
the  high  energy  tail  of  the  electron  distribution  is  enhanced  by  the 
electric  field,  leading  also  to  runaway  electrons.  These  calculations, 
however,  ignore  inelastic  collisions  and  the  electric  field  effects  on 
ionization  balance.  Calculated  runaway  electron  production  rates  also 
agree  poorly  with  experiment.  Finally,  recent  work  on  the  effects  of  steep 
gradients  on  heat  flow  have  shown  that  the  process  is  nonlocal,  i.e.,  no 
longer  described  by  the  local  gradient  in  temperature,  and  that  the  high 
energy  tail  of  the  electron  distribuiton  may  either  be  depleted  or 
populated  relative  to  a  Maxwellian  by  this  energy  flow. 

III.  APPROACH 

As  was  done  in  the  determination  of  the  ionization  balance  and 
radiation  loss  rates  of  different  atomic  weight  plasmas,  the  problem  of 
calculating  the  electron  energy  distribution  function  is  being  approached 
as  one  involving  moderate-Z,  partially  ionized,  non-LTE  plasmas.  Thus,  a 
model  is  being  developed  that  will  allow  the  Z-scaling  b^avior  of  plasma 
conductivities  and  plasma  ionization  balance  to  be  calculated  as  a  function 
of  the  strong  and  rapidly  changing  S  and  3  fields  of  an  imploding  Z-pinch. 


To  begin  the  electron  calculation,  atomic  models  of  hydrogen  and  argon 
plasmas  are  being  developed.  The  hydrogen  model  is  important  for  several 
reasons.  One,  it  is  an  important  dilutant  to  be  used  in  the  study  of  the 
implosion  dynamics  of  mixed  element  loads.  It  is  ionized  quickly  during 
the  initial  phases  of  the  implosion  and  provides  an  additional  sources  of 
electrons  for  plasma  heating  and  for  enhancing  excitation  and  ionization 
rates  of  the  moderate-Z  elements  in  the  mixture.  Two,  a  hydrogen  model  is 
easily  scaled  for  use  in  the  development  of  K -shell  models  for  moderate-Z 
plasmas.  And  three,  it  is  an  atomic  system  for  which  much  is  known  and  for 
which  conductivity  and  ionization  balance  calculations  can  be 
benchmarked.  The  argon  model,  on  the  other  hand,  will  be  developed  to 
demonstrate  that  the  calculation  will  also  work  in  general  for  moderate-Z 
elements  and  to  determine  Z-scaling  effects  that  will  occur  because  of  the 
different  non-LTE  behaviors  of  different  ionization  stages.  The  excitation 
of  the  argon  system  will  also  show  different  transient  ionization  and 
conductivity  behavior. 

The  calculation  of  electron  distributions  from  inelastic  scattering 
processes  entails  a  new  approach  to  the  problem  of  ionization  equilibria  as 
well.  One  must  work  with  cross-sections  rather  than  rate  coefficients, 
which  are  to  be  obtained  in  this  case  as  cross-section  averages  using  the 
actual  rather  than  an  assumed  Maxwellian  electron  distribution.  Moreover, 
the  cross-sections  must  be  expanded  in  the  same  way  as  f,  i.e.,  their 
angular  distribution  must  be  known  in  order  to  carry  out  the  conductivity 
calculation. 

IV.  FUTURE 

The  main  goal  of  the  present  effort  is  to  self -consistently 
incorporate  the  effects  of  inelastic  collisions  in  the  calculation  of  the 
shape  of  the  distribution  function  f^  and  on  the  ionization  balance  and 
radiation  loss  rates.  A  conductivity  calculation  will  then  be  set  up  at 
first  using  this  non-Maxwell i an  fg  and  then  setting  up  a  calculation  that 
includes  the  simultaneous  effects  of  a  strong  E  field  and  inelastic 
collisions  on  fg  and  f-, .  These  calculations  will  be  applied  to  an 
investigation  of  the  power  flow  problem  in  imploding  z-pinchs.  These 
calculations  will  be  applied  mainly  to  critique  implosion  behavior  in 
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